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Capítulo 1 
Movimiento ondulatorio 


1.1. Concepto de onda que se propaga. Frente de onda. 


Se denomina onda a cualquier perturbación del equilibrio que se propaga de un punto a otro del espacio. Las ondas están 
presentes en multitud de fenómenos de la vida cotidiana. Podemos, en una cuerda tensa retorcer un extremo o dar un tirón, 
o producir una variación brusca de presión en un punto de un gas, y esta perturbación se propagará a lo largo del sistema. 
Por ejemplo, el sonido de la voz humana no es más que una perturbación de las condiciones de equilibrio producida en las 
cuerdas vocales, que se propaga por el espacio y es detectada en el oído de otra persona. Las ondas luminosas también son 
una perturbación de las condiciones de equilibrio que se propaga por el espacio. Otro ejemplo de ondas serían las olas de 
agua que se acercan a la costa o las pequeñas ondas que se propagan por la superficie del agua. Nótese que en ningún caso 
de las ondas mecánicas que hemos visto la materia viaja, sino solo su estado de movimiento, hay transporte de energía y 
momento en el espacio, pero no de materia. 


Se denomina frente de onda al lugar geométrico de los puntos del espacio sometidos al mismo estado de perturbación 
(La perturbación en todos esos puntos es la misma y tiene la misma fase). La dirección de propagación siempre es 
perpendicular al frente de ondas. Se denomina a la onda onda plana cuando el frente de ondas es plano. Dependiendo 
de la forma del frente de onda las ondas pueden también ser esféricas, cilíndricas, etc ... Las ondas planas y las esféricas 
tienen una especial importancia y por ello se les dedicará un mayor espacio. 


Asimismo distinguiremos entre ondas longitudinales y transversales. En las ondas longitudinales, la dirección en la cual 
varía la magnitud que define la perturbación coincide con la dirección de propagación de la misma, mientras que en 
las ondas transversales, la dirección de variación de esta magnitud es normal a la dirección de propagación. El sonido 
propagándose en un gas es un ejemplo de onda longitudinal, mientras que una onda electromagnética o una onda que se 
propaga a lo largo de una cuerda tensa son ejemplos de ondas transversales. En algunos casos estos dos tipos de ondas 
van combinados, como ocurre con las ondas sísmicas. 


1.2, Descripción matemática de la propagación para una onda plana 


Consideremos una función representada por la curva central de la figura 1-1. Reemplazando la variable x por x — a con 
a una cantidad positiva, se tiene la función de la derecha. La forma de la función no ha cambiado, pero sí dónde está su 
origen. Se tienen los mismos valores de la función, pero para puntos desplazados hacia la derecha una cantidad a. De 
forma similar se tiene que al reemplazar en la función original x por (x + a), la función es la misma, pero desplazada 
hacia la izquierda una cantidad a. 


fa+a) 1) fa-a) 
- Ad : Mm 
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Figura 1-1 Descripción matemática de la propagación de una onda 


Si en vez de la constante positiva a hubiéramos utilizado ct, donde t es el tiempo y e una velocidad, se tendría que f(x—ct) 
correspondería a la curva original, pero viajando sin deformación hacia la derecha a una velocidad c. Análogamente, 
F(z + ct) correspondería a la función original desplazándose sin deformación hacia la izquierda con una velocidad —c. 
Denominamos perfil de onda a la representación de f(x + ct) en un instante dado. 
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Así, una expresión matemática de la forma 

x= Fx + ct) (1-1) 
es apropiada para describir una perturbación que se propaga sin deformación en el eje x, en los sentidos positivo [ f(x — 
et) ] y/o negativo [ f(x+ct) ]. Este tipo de movimiento se denomina movimiento ondulatorio. A las ondas que se propagan 
se les denomina también ondas progresivas en contraposición a las ondas estacionarias que se estudiarán más adelante. 


1.3, La ecuación de ondas para las ondas planas en una dimensión 


Ya sabemos la forma general que debe tener una onda que se propaga en una dimensión. Sin embargo, a la hora de saber 
como se propaga una perturbación y concretamente con qué velocidad, la expresión (1—1) anterior no es práctica ya que 
es difícil de obtener a partir de las características físicas del proceso. Afortunadamente es fácil obtener una relación que 
debe cumplir la solución (1-1), y que está más relacionada con las características físicas del proceso de propagación. Para 
ello utilizaremos la variable u definida como el argumento de x, es decir u = x + ct y derivaremos la ecuación anterior 
respecto del tiempo y de la posición a través de esta variable u, quedando 


Ox _ 0x0 Ox xd 0 
0 — dudt du 02" 90 du? 


Ox 0x0u 0x Px  0xd0u _0x 


a A CE AAA A 
0x  0udz du 012 due%tdx du 
que igualando las segundas derivadas parciales de x respecto de u da inmediatamente la expresión 


Px _ ¿Ox 
ot? Ox? 
conocida como ecuación de ondas o ecuación de D”Alembert, que debe cumplir toda onda plana que se propaga con 
velocidad ce y sin deformación en una dimensión, y cuyas soluciones deben ser obviamente de la forma x(x,t) = x(u+ct). 


Nótese que aunque hemos obtenido la ecuación de D” Alembert a partir de la función x(x + ct) se obtendría el mismo 
resultado para la función x(x — ct). 


(1-2) 


La ecuación de D'Alembert (1-2) es una ecuación lineal en x, ya que ni x ni sus derivadas aparecen más que como 
expresiones de primer grado. Así, se cumple el principio de superposición y cuando a un punto llega simultáneamente un 
conjunto de perturbaciones el resultado de dichas perturbaciones en dicho punto se puede obtener como la superposición 
de los resultados producidos por cada una de las perturbaciones individualmente, y por tanto la función que describe la 
perturbación completa como la suma de funciones que describen cada una de las perturbaciones individuales. 


En la siguiente sección veremos la gran utilidad de esta expresión, obteniendo las velocidades de propagación de distintas 
perturbaciones en distintos medios. 


1.3.1. Análisis de Fourier de una onda periódica. 


Fourier demostró que las ondas periódicas de formas complejas pueden representarse como suma de una serie de ondas 
armónicas (ondas representadas por funciones seno o coseno). Consideremos una función x(t) que represente la magnitud 
que tiene en una cierta posición una onda periódica de periodo 7". Según el teorema de Fourier esta función puede 
representarse mediante una serie de la forma 
x(t) = ay + ay coswt + b¡ senwt+ (1-3) 
+ 49 C08 2wt + ba sen 2wt-+ 
+ ... + 
+ 4, cos nut + b,, sen nuvt 


donde w = 27 /T. 


Dependiendo del perfil de la onda, el número de términos necesario para obtener una buena representación de la onda 
original será mayor o menor. Se denomina frecuencia angular fundamental o armónico angular principal a la frecuen- 
cia angular w y a las frecuencias angulares 24, 3w,...,nw armónicos angulares superiores o sobretonos angulares. Se 
denomina análisis de Fourier al proceso de determinación matemática de los coeficientes a; y b;. 
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Figura 1-2 Análisis de Fourier de una onda cuadrada, conteniendo los tres primeros términos no nulos del desarrollo en serie de Fourier. Superpuesto 
a la onda cuadrada se muestra la suma de estos tres primeros términos del desarrollo. 


Es decir cualquier movimiento periódico puede analizarse como una superposición de funciones armónicas. Por ejemplo, 
en la figura 1-2 se muestra el caso concreto de una onda cuadrada con los primeros términos no nulos de su descomposi- 
ción de Fourier (Obsérvese que los términos pares deben ser nulos). 


No es posible analizar un movimiento no periódico mediante una serie de Fourier. Sin embargo sí es posible hacerlo 
utilizando una integral de Fourier. 


1.3.2, Propagación de una onda armónica plana 


Como se acaba de ver al tratar el análisis de Fourier, un caso de especial importancia se tiene cuando la forma de la 
función (1-1) es armónica de la forma 


x(x,t) = xosen k(x + ct) (1-4) 


donde k no es simplemente una constante que hace adimensional el argumento del seno, sino que tiene un significado 
especial. Consideremos la ecuación anterior, pero reemplazando z por x + 27/k, 


2 2 2 
X ( + q) =N (=+ 7 ta) = xosenk (2. ta) = xosen[k(xw + ct) + 2] = x(u + ct) 


Se ve claramente que 
A=21/k (1-5) 


es el periodo espacial de la onda dada por la función (1-4), es decir, la función se repite a sí misma cada longitud A, 
denominándose a esta cantidad A longitud de onda. k = 21 /A representa el número de longitudes de onda en la distancia 
2r y recibe el nombre de número de onda angular , aunque dada su importancia en los fenómenos de propagación nos 
referiremos a él como número de onda.? Sus unidades son radianes de fase por unidad de longitud (o ciclos en cada 271 


A veces, y especialmente en textos de espectroscopía, se denomina número de onda a la magnitud k = 1/A = k/21r, también denominado 
número de onda espectroscópico , que corresponde al número de longitudes de onda por unidad de longitud (número de ciclos por unidad de longitud), 
y que tiene unidades de inversa de longitud. Este es análogo, para oscilaciones en el espacio, a la frecuencia y para las oscilaciones temporales. 
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unidades de longitud) y es análogo, para la periodicidad espacial, a la frecuencia angular w en la periodicidad temporal, 
representando la tasa de cambio de la fase total con la distancia. 


Así, podemos representar una onda plana sinusoidal que se propaga hacia la derecha a lo largo del eje x mediante una 
ecuación de la forma 


x(x,t) = xo sen k(x — ct) = xo sen lo — ct) (1-6) 
Esta ecuación puede también escribirse en la forma 
x(x,t) = xosen(kz — wt) (17) 
donde mr Pre le 
A 


es la frecuencia angular de la onda. Como la frecuencia angular está relacionada con la frecuencia y con la que la situación 
física varía en un punto, a través de la expresión w = 27p, se tiene la relación 


277 w 
as de (1-9) 


entre la longitud de onda A, la frecuencia y y la velocidad de propagación de la onda c. 


Podríamos haber escrito lo anterior en función del periodo Tper = 277 /w = 1/v en la forma 


x(x,t) = xo sen 277 G E ) (1-10) 


Como se ve en las ecuaciones anteriores, para una onda armónica plana que se propaga se tiene una doble periodicidad, 
espacial y temporal. En un instante de tiempo dado, como si hubiéramos hecho una fotografía instantánea de la pertur- 
bación, hay una periodicidad espacial en la onda, ya que ésta se repite cada longitud de onda A. Igualmente, en un punto 
fijo del espacio afectado por la perturbación al observar la evolución de la perturbación en ese punto se ve que hay una 
periodicidad temporal, ya que la onda se repite cada periodo. 


1.4. Propagación de ondas transversales en una cuerda tensa. 


1.4.1. Ecuación de onda y velocidad de propagación 


Consideremos ahora una cuerda que se encuentra tensa bajo la acción de una fuerza tensora F. En el equilibrio mecánico 
la cuerda se encuentra prácticamente recta (salvo una pequeña catenaria que no consideraremos). Si en un punto de la 
cuerda desplazamos ésta perpendicularmente, un segmento dx que se encuentra en un punto zx se habrá desplazado en 
la dirección perpendicular una cantidad u, como se ve en la figura 1-3, manteniéndose la fuerza tensora FF la misma 
en ambos extremos. Esta fuerza es debida a la tracción que realiza el resto de la cuerda a ambos lados del segmento, 
pero debido a la curvatura de la cuerda ambas fuerzas ya no se aplican en la misma dirección y sentidos opuestos. Las 
componentes verticales de cada una de las fuerzas son 


F,y = —Fseno 
s F, = Fseno!, 


por lo que la fuerza resultante sobre el segmento AB de la cuerda es 
F, = F(sena” — sena) 
Como la curvatura es en general pequeña es posible aproximar los senos por las tangentes correspondientes, quedando 


Otana 


Ox 


tan a es la pendiente de la curva formada por la cuerda en el punto x, y puede expresarse como 9u/0.x, quedando para la 
fuerza neta en la dirección y 


F,=F(tano—tana) —> dF,=YFdtena=J dx 


dE, =F—- dx = F- dr. (1-11) 
x x 
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Figura 1-3 Sistema de fuerzas que actúa sobre un elemento de cuerda de longitud dx que se ha desplazado en la dirección perpendicular una distancia 
u. 


Esta fuerza neta debe ser igual al producto de la masa del segmento de cuerda, px, por la componente de la aceleración 
del elemento de cuerda en la dirección y, 
u Pu 
donde p, es la densidad lineal de masa de la cuerda. De las dos ecuaciones anteriores se tiene 
Pu FO ¿Ou 


92 por "da 


que tiene la forma de la ecuación de ondas, lo que quiere decir que una perturbación transversal en la cuerda se propaga 
longitudinalmente a lo largo de la misma con una velocidad 


(1-13) 


PR (1-14) 


La fuerza resultante paralela al eje x es 
F, =F cosa! — Feosa = F(cosa! — cos a) 


que cuando el ángulo es muy pequeño es, con muy buena aproximación, nula. Nótese que en los cálculos anteriores se ha 
considerado únicamente el caso en el que la amplitud es pequeña. 


1.4,2, Energía de una cuerda vibrante 


Calculemos ahora la energía de una cuerda vibrante. La energía cinética de cada elemento de la cuerda de longitud original 
dx como el que vimos en la figura 1-3 vendrá dada por 


1 gu? 1 du” 


donde u es el desplazamiento de la cuerda y p, su densidad lineal de masa. Por otra parte, dicho elemento de la cuerda es 
estirado hasta una longitud dx” dada por 
E E 
2 1 0x 


2 2 
de! = y) de? + (5) da? =dx4/1+ (5) = dx 


por lo que cambia su longitud en 
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bajo la acción de la fuerza tensora F. El trabajo realizado nos dará la energía potencial del elemento, 


2 
dE, = F(dx' — dx) = Z (5) dz. (1-16) 


La energía total por unidad de longitud de la cuerda será 


dE 1 (0Y 1 2¿(0Y'_1 2|/(0uY”, (10uY 
de = 30 (5%) tala) 2 ll TA y 
F 


y la energía total será la integral de esta expresión extendida a toda la cuerda. 


Densidad lineal de energía para una onda armónica 


Si por la cuerda se propaga una onda armónica de la forma u = uy cos(kx — wt) la expresión anterior da 
dEsot 
dx 


y su promedio temporal nos da la densidad lineal de energía promedio 


= pic“uzk? sen? (kx — wt) 


lo222_1 232 
Et = gane uyk? = ¿Pu (1-18) 
para una cuerda vibrante a través de la que se propaga una onda armónica. 
Potencia transmitida por una onda armónica a través de la cuerda. 
Consideremos el extremo izquierdo del elemento dx de la figura 1-3 como si fuera el extremo izquierdo de la cuerda, que 


se está desplazando con una velocidad ú bajo la acción de la componente vertical de F. La potencia transmitida a través 
de la cuerda desde ese extremo vendrá dada por 


que en función de la fuerza tensora F queda 
O o 
dt Ox dt dx” 


donde el signo — aparece como consecuencia de que si la pendiente du/0x es positiva la fuerza recuperadora va hacia 
abajo, y de que en una cuerda tensa F es siempre positivo. 


Si la onda que se propaga por la cuerda es armónica de la forma u = uy cos(kz — wt), ésto da 
dW , 
a * Fkwuz sen*(kx — wt) = cpjw*uj sen*(kz — wt) 


con F = pic? y k = w/c. Promediando temporalmente esta expresión y haciendo uso de la ecuación (118) se obtiene la 
potencia promedio que se propaga a través de la cuerda con una onda armónica, 


dW 1 1 
(d)> cpu” = el gouotaz) = Ep. (1-19) 


1.5, Propagación de ondas transversales en una membrana tensa 


Analizaremos ahora la propagación de ondas transversales en una membrana tensa, como la superficie de un tambor o 
el diafragma de un micrófono. En nuestro análisis supondremos que la membrana es delgada y uniforme y que tiene un 
comportamiento perfectamente elástico, con una rigidez muy baja y sin disipación de energía por rozamiento, siendo ps 
la densidad superficial de masa de la membrana. 


Nos restringiremos a desplazamientos pequeños y consideraremos que la fuerza tensora por unidad de longitud F, es 
uniforme en la membrana e isótropa. Así, la fuerza ejercida a través de un segmento de longitud dl por el material de un 
lado sobre el del otro es F¡ dl. 


1.5. Propagación de ondas transversales en una membrana tensa L 


1.5.1. La ecuación de onda para una membrana tensa 


Utilizaremos coordenadas cartesianas para resolver el problema general. Nótese que el uso de un sistema de coordenadas 
adecuado en un problema particular puede simplificar enormemente el tratamiento del mismo. Si denominamos zy al 
plano de la membrana y z al eje en el que tienen lugar los desplazamientos (a los que denominaremos u), la fuerza 
recuperadora que actúa sobre un elemento de superficie infinitesimal dS = dx dy será la suma de las componentes según 
z de las fuerzas recuperadoras en los lados paralelos a x e y. Mediante un razonamiento análogo al seguido para la cuerda 
tensa en la sección 1.4, la fuerza neta en la dirección z que actúa sobre el elemento de superficie d,S debido a las fuerzas 
en la dirección y es 
O tan Qy 0 du q Hu 
—_d dy — Tdzx = 

Ox Ox Ox Fog a 
mientras que la fuerza neta en la dirección z que actúa sobre el elemento de superficie d.S debido a las fuerzas en la 
dirección x es 


Fi dy = Fi — dedy 


q Hu 
> ó A dedy. 


La suma de ambas fuerzas debe ser igual al producto de la masa del elemento de superficie por su aceleración en la 
dirección z, 
gx? * Oy? 


que nos da la ecuación de ondas para una membrana elástica tensa, 


2 
(+ ex e a) dady = pudady 2 (1-20) 


2 2 2 2 
9 un q 200_ 1050 (1-21) 


donde 


N á (1-22) 


es la velocidad de propagación de las mismas. El primer miembro de la ecuación (1-21) es la laplaciana de u en coor- 
denadas cartesianas planas, y en esa forma es muy útil para analizar membranas rectangulares. La ecuación de ondas así 
expresada es válida para cualquier sistema de referencia en el que los ejes están rotados respecto de los anteriores en el 
plano xy. 


Si la membrana es circular como el parche de un tambor, es conveniente expresar la laplaciana en coordenadas polares 
planas (r, 0), con  = r cosd y y = r sen Ó, quedando la ecuación de ondas (1-21) de la membrana en la forma 


du 10 10%u 10 
Hr poro Ss a e pr 23 
Or? " rór  r20902 e 9t (4-23) 


1.5.2. Solución de ondas planas en coordenadas cartesianas 


Buscaremos ahora una solución de la ecuación (1-21). Para ello supondremos que la solución general puede expresarse, 
utilizando el método de separación de variables, en la forma 


u(z, y,t) = X(1)Y (yT(t) (1-24) 


donde X (x=) es función exclusivamente de x, Y (1) es función exclusivamente de y y T'(t) es función exclusivamente de 
t. Sustituyendo en la ecuación (1-21), se tiene 


2 2 
A A 
que, dividiendo por u(x,y,t) = X(1)Y (y)T (t) queda 
1 d2X(x) 1 4Y(yY_1 1 d7(t_ _p2 (1-25) 


X(x) di? “Y(y) d?  ¿2T(t) de 
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donde, al ser los dos miembros de la ecuación función de distintas variables independientes, y al deber mantenerse la 
igualdad para cualquier valor de dichas variables, ambos deben iguales a una constante a la que denominaremos —k?. Así 
tendremos por una parte la ecuación diferencial de la dependencia temporal, 


d?T (t) GTO. : 
2 mo 2 — “e 1 - 
e? T(t) dt? Ñ di? (kc)*T (t) de w"T(t)=0 (1-26) 


con w = kc, que es una ecuación similar en forma a la del oscilador libre que se ve en la sección 2.2.1 de la referencia [4] 
y cuya solución es 


T(t) = 444%) (127) 


y por otra parte la ecuación diferencial de la dependencia espacial, 


1 d?X(x) 1 YY _ , 
X(0) da? =P Y) dy? = —k (1-28) 


Cada uno de los términos del primer miembro de la ecuación depende de una única variable independiente y su suma debe 
ser constante. Como ésto debe ser así para cualquier valor de las variables x e y, ambos deben ser constantes, por lo que 
si denominamos —k? y E a las constantes, se tiene que las correspondientes ecuaciones diferenciales separadas, junto 
con sus soluciones, son 


1 d2X(x) hera) 
- = = ds z 1-2 
Xx de +k2X(2) =0 > X(12) = A, € (129) 
1 dy (y) 2 (k,yY—0y) 
ds 28 Y Y 3 
YO ar + k2Y (y) =0 > Y (y) = Ay € (1-30) 


ya que cada una de las ecuaciones diferenciales es de nuevo de la forma de la ecuación diferencial del oscilador libre. 
Asimismo se cumple 


pep poÉ (131) 
z y e p 


La solución completa para el desplazamiento en la dirección z queda entonces 
ula y t) . Ae itr+tky=wt-a) ES pen Er-ut-a) (132) 


con == (1, y) y k= (ki, ky), donde hemos puesto los signos así por comodidad, ya que la misma ecuación con +k, Ek, 
y +w es también solución de la ecuación (1-21). Esta solución corresponderá a una onda que se propaga sin deformación 
con velocidad c = w/k en la dirección dada por k. 


Analicemos un poco más la solución obtenida. Para ello observemos el frente de onda, lugar geométrico de los puntos en 
los que la perturbación tiene el mismo valor, es decir en los que la fase total es constante e igual a kx + kyy — wt— o. 
Como a es constante, kx + kyy — wt también será constante en' el frente de onda e igual a un cierto valor al que 
denominaremos ¿, que será característico de cada frente de onda. En un instante genérico t los puntos del frente de onda 
satisfarán la ecuación 


kg2 + kyy —wt = 0, (1-33) 


es decir, en cada instante el frente de onda es una línea recta en el plano xy de la membrana, dada por y = ($ + wt)/ky — 
ky3/k,, que es perpendicular a k. Así, la ecuación (1-32) representa una onda plana bidimensional que se propaga en la 
membrana en la dirección k. 


Si en un mismo instante consideramos el frente de onda anterior y el siguiente, correspondiente a un aumento de 27 en la 
fase, se tendrá, con 7 y Fa, 


7 =0+wt OT 
> k-(m—F,)=2=Ak 1-34 
E-Fo=ó+wt+2m si si 


1 


donde A es la longitud de onda, que se obtiene de la diferencia de las proyecciones de 7 y Fa sobre la dirección de 
propagación y que muestra la distancia entre dos frentes de onda consecutivos y k = ,/k2? + k$ = w/c es el número de 


onda de la onda plana, mientras que a k se le denomina vector de onda. 
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1.5.3. Solución en coordenadas polares planas 


Si resolvemos el problema anterior en coordenadas polares planas con la ecuación de ondas (1-23) y utilizamos el método 
de separación de variables probando con una solución del tipo 


u(r,0,t) = R(rjO(0)T (t) (1-35) 
obtenemos para la parte temporal una dependencia similar a la que obteníamos en coordenadas cartesianas, 
Tb) 4 EUA (1-36) 
quedándonos el resto como la denominada ecuación de Helmholtz en coordenadas polares 


p2 .s o) + hi2 L_40(0) (1-37) 


R(TA dr? dr 0(0) de? 


donde k = w/c. El primer miembro de la ecuación es una función únicamente de r mientras que el segundo miembro es 
una función únicamente de 0 y ésto debe ser así para cualquier valor de r y 9. Esto solo es posible si ambos miembros son 
constantes. Si denominamos m? a esta constante el segundo miembro queda como 


d20(0) 2 
——_—— == A) 
49 +m*0(0) =0 (1-38) 
que tiene una solución armónica de la forma 
O(0) = Ap cos(mb — yop). (1-39) 


Para valores de r y t dados, esta función debe repetirse cada 27 radianes. Para que ello sea posible, m no puede tener 
valores cualesquiera, sino únicamente valores enteros positivos o 0, m = 0, 1,2, 3... de forma que se cumpla siempre que 
0(0) = O(0 + 27). Con estos valores posibles de m se obtiene, a partir del primer miembro, la ecuación de Bessel 


PR 1dR 3 
Ea e 2 + (1 de =) R(r) =0 (1-40) 


Como es una ecuación de segundo orden debe tener dos soluciones linealmente independientes para cada valor del pará- 
metro m. Estas soluciones son dos funciones especiales llamadas funciones de Bessel de primera especie Jm(kr) y 
funciones de Bessel de segunda especie Y, (kr) ambas de orden m. 


Las funciones de Bessel de segunda especie también son denominadas funciones de Neumann y se hacen infinitas para 
kr = 0 aunque tienen valores finitos en otro caso. Por ese motivo no pueden representar un posible desplazamiento 
de una membrana circular completa, que es el caso más habitual en el que se usan coordenadas polares, por lo que en 
muchos casos bastará con utilizar las funciones de Bessel de primera especie. Sin embargo, puede ser necesario utilizarlas 
en otro tipo de problemas, como en el caso de una membrana anular (para la que obviamente no hay problemas en su 
centro). Todas las funciones de Bessel de primera especie se hacen O en el origen salvo Jo que se hace 1 y oscilan con 
amplitudes que disminuyen según aumenta kr de una forma que recuerda a la de un oscilador amortiguado. Las funciones 
de Bessel están tabuladas en los manuales habituales de fórmulas matemáticas y definidas en muchos programas de cálculo 
matemático, aunque con frecuencia solo las dos primeras de cada tipo. Ésto es debido a que son recurrentes y basta con 
conocer dos sucesivas para conocer todas las demás utilizando la fórmula de recurrencia 


1.6. Propagación de ondas mecánicas en sólidos. 


En los sólidos se tienen fuerzas recuperadoras elásticas frente a deformaciones tanto longitudinales (de tracción y compre- 
sión) como transversales (cortantes). Esto hace que en un sólido puedan aparecer tanto ondas longitudinales como trans- 
versales. Dentro de éstas últimas se encuentran las ondas de torsión, producidas por fuerzas recuperadoras de torsión. 
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0 10 í 


Figura 1-4 Representación gráfica de las tres primeras funciones de Bessel de primera especie (arriba) y de segunda especie (abajo) 


1.6.1. Propagación de una onda mecánica en un sólido elástico isótropo e ilimitado. 


Para obtener las características de la propagación de una onda en un sólido elástico isótropo e ilimitado es necesario 
partir de la ecuación de movimiento de un elemento de volumen infinitesimal del sólido, en la que se tienen en cuenta las 
acciones de las fuerzas elásticas a través de las superficies del elemento, 


dE, _ 00%x _ OTzy E OTaz A Ñ zx 
E TN E E 
dFy _ Otys , 00yy , OTyz - Pr 
e Td de e da 


dE, _ Ox , OToy , D0%% E oz 
dV dr "y "o PR Pa 


donde y = (9., 9y, 92) representa la fuerza por unidad de masa debida a la acción de fuerzas de masa como la gravitatoria 
y las otras componentes representan la acción de fuerzas elásticas. La expresión anterior puede escribirse de una forma 
más compacta en la forma siguiente, en la que además hemos despreciado el efecto de la fuerza gravitatoria, 


Pz =V6 (141) 


donde d es el tensor de tensiones y Y el vector de desplazamientos, siendo 7 = Fy + U, con Fo la posición del elemento en 
el equilibrio. 


Consideremos la deformación del elemento prismático infinitesimal de dimensiones dx, dy y dz a lo largo de los ejes 
coordenados. Como consecuencia de la perturbación estas dimensiones pasan a ser du + du, dy + du, y dz + du,, donde 
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Uy, Uy y u, SON los desplazamientos en cada una de las direcciones de los ejes coordenados y, al estar las aristas dirigidas 
originalmente según los ejes, 


JU du 0 
dus = te, duy = dy . Uly a nda 
por lo que la deformación cúbica o deformación unitaria de volumen queda 
dv 0 0 0) > 
NA E MR (1-42) 


es decir, 0 es igual a la divergencia del vector desplazamiento en ese punto. 


Las relaciones entre tensiones y deformaciones vienen dadas por la conocida ley de Hooke 


vE Es E Es 
= Vu I + ——é = AL Vul + 2Gé = - E 143 
ó a+ Ha + an ¡Vul + 2Gé = 101 + 2Gé ( ) 


donde » es el coeficiente de Poisson, Y el módulo de elasticidad longitudinal o módulo de Young y é el tensor de defor- 


maciones a a F 
0Ux OU QUe 


dx 0y 0z 
a 0%, 0uUy 0QUuy 
¿== | == e 1-44 
ú 0x 0y 0z |” ( ) 
dia dia Pte 
dx 0y 0z 
mientras que A; y G son las constantes de Lamé, 
vE 
M = 3 145 
'= 114001 - 2) si 
s E 
A 1-45b 
E 2(1 +.) ) 
denominándose a G' módulo de rigidez. Combinando las ecuaciones (1-41) y (1-43) se obtiene 
0% vE e 2 — HP 0 
=- DY + VW UU =AV(Vu) + 2G Vu 1-46 
Poe A+DNad— rl úl (1 +») 7 ida e qe 


Si las deformaciones producidas por la onda son pequeñas, las fuerzas recuperadoras que intervienen pueden considerarse 
elásticas. En esta expresión hemos utilizado la notación 


<<: 
mm 
li 


o (1-47) 


Estudiaremos en primer lugar una onda plana que se propaga en un sólido isótropo ilimitado!, es decir una onda en la 
que los desplazamientos ú son función únicamente de una de las coordenadas, por ejemplo zx, y del tiempo. Así todas las 
derivadas parciales respecto de y y 2 desaparecen y la ecuación anterior se puede, para los distintos ejes, descomponer en 
las ecuaciones 


Pus ¿Uy 

az * Ci 792 (148a) 
uy ¿Ot 4 
gr = y 5) 2 (1-48b) 
du, O n 
qa aa” (1480) 


testo es importante, porque en estas condiciones apenas hay deformación transversal para una longitudinal dada 
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Las ecuaciones (1-48) tienen la forma de la ecuación de onda para ondas planas (1--2) que se propagan sin deformación 
en la dirección x con velocidades c, o c+. La primera de las ecuaciones correspondería a una onda longitudinal (produce 
un desplazamiento a lo largo de la dirección de propagación x) que se propaga con velocidad 


Ev 2G_ E( —v) 


E CIN TN AREA ON 1-40 
p(1+vw)(1 — 21) + P p(+»w)(1 — 21) 1 


q 
donde G es el módulo de rigidez (1-45b) y p la densidad del medio. Las otras dos ecuaciones corresponderían a ondas 
transversales (producen desplazamientos en las direcciones y o z perpendiculares a la dirección de propagación) que se 
propagan con velocidad 

E G 


E A -5 
AE ae 


> 


De esta forma vemos que una onda elástica general que se propaga en un sólido isótropo ilimitado en realidad consta de 
dos ondas que se propagan independientemente una de otra, una de ellas longitudinal (con velocidad c,) y otra transversal 
(con velocidad c;). La velocidad de propagación es distinta según la perturbación causada por la onda sea longitudinal 
o transversal. En los sólidos isótropos, la velocidad de las ondas longitudinales es mayor que la velocidad de las ondas 
transversales (ya que, normalmente, 0 < y < 1/2), siendo la relación entre ambas 


o O TA 
1-v(1421) Y 1-2" 


(1-51) 


Un sismólogo puede estimar la distancia al epicentro de un temblor midiendo la diferencia en los tiempos de las ondas 
longitudinales y transversales. 


En los sólidos anisótropos, las propiedades elásticas no son las mismas en todas las direcciones. Por tanto, la velocidad de 
las ondas longitudinales y transversales depende de la dirección de propagación y en el caso de las ondas transversales, 
también de su polarización, esto es, de la orientación del plano formado por la dirección de propagación y la de variación 
de la perturbación, al que se conoce como plano de oscilación. 


1.6.2. Propagación de ondas elásticas longitudinales planas en una barra fina. 


Las ondas que se propagan en una lámina delgada o en una barra fina son a menudo distintas de las que se propagan en un 
sólido ilimitado. Esto es así cuando la longitud de onda de las ondas es grande comparado con el espesor de la lámina o el 
diámetro de la barra. En el caso contrario, en el que las longitudes de onda son pequeñas frente al espesor del elemento, 
el medio puede considerarse en la práctica como ilimitado en todas las direcciones y las relaciones (1-48) siguen siendo 
válidas. 


in el caso de una barra fina, le deformación que se produce en la dirección longitudinal no depende de ninguna fuerza 

de tracción o compresión transversal y por lo tanto se trata en la práctica de una tracción o compresión simple, donde la 
ley de Hooke es válida en su forma más sencilla, v = He. Si no se tratara de una barra fina, aparecerían también fuerzas 
recuperadoras elásticas de tracción o compresión en las direcciones perpendiculares a la de propagación y deberemos 
utilizar las ecuaciones de la sección anterior. El caso de la lámina delgada es un caso intermedio entre ambos. 


Obtención de la ecuación de ondas 


Consideremos una barra (fina) y perturbemos uno de sus extremos, por ejemplo golpeándola. Esta perturbación se propaga 
alo largo de la barra y finalmente se siente al otro extremo. La onda elástica producida por la perturbación se ha propagado 
a lo largo de la barra. Estudiemos esta propagación. Para ello consideraremos una barra de sección transversal uniforme, 
sometida en caras opuestas de cada una de sus secciones a una fuerza F"(1) que en general será distinta a lo largo de la 
barra. 


Consideremos ahora las fuerzas que actúan sobre una rodaja de un cierto espesor. En una rodaja como la que se muestra 
en la figura la magnitud de la fuerza a un lado y a otro tendrá en general distinto valor al ser F(x) distinta a lo largo 
de la barra. Sin embargo, si consideramos la fuerza que se ejerce por la rodaja sobre el resto del material en una de las 
bases de la misma y la que ejerce el resto del material sobre la rodaja a través de la misma cara, por el principio de 
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dx 
SS E A 


T+U 


Figura 1-5 Propagación de una perturbación a lo largo de una barra fina. Una rodaja de espesor dx situada en x se desplaza a 2/ = u + 2 durante el 
paso de la onda y su espesor pasa a ser da” = dz + du. A través de cada una de sus caras actúan unas fuerzas FF y F”. 


acción y reacción ambas tendrán la misma magnitud, pero distintos sentidos. Para describir como es la fuerza de tracción 
o compresión que actúa en la barra en una determinada sección utilizaremos este conjunto de fuerzas a través de una 
sección al que denominaremos .F. El valor de F en esa sección será |F(x)| a través de la sección y su signo positivo si 
produce una tracción y negativo si produce una compresión, y en el caso general variará de un punto a otro de la barra. 


Bajo la acción de estas fuerzas, una sección genérica de la barra que originalmente se encuentra en una posición x sufre un 
desplazamiento u paralelo al eje, pasando a la posición x + u. Si este desplazamiento es el mismo en todos los puntos de 
la barra no se produce deformación alguna, sino un desplazamiento de toda la barra como un sólido rígido. Consideremos 
ahora el caso en el que se produce deformación. Para ello estudiaremos qué ocurre con la rodaja de espesor di que 
originalmente está en x. Ya hemos visto que la cara izquierda se desplaza a 1 + u. La cara de la derecha se desplazará en 
u! = u + du, pasando su espesor de du a de” = dx + du. El alargamiento (o acortamiento si du es negativo) de la rodaja 
ha sido por tanto du. 


Utilizaremos la deformación unitaria (deformación por unidad de longitud) en la dirección x, definida por 


e= se (1-52) 
y la tensión elástica, definida por 
o=2 (1-53) 


sabiendo que en un estado de tensiones unidimensional y pequeñas deformaciones, ambas están relacionadas por la ley 
de Hooke correspondiente a una tracción o compresión simple 


a =E€ (1-54) 
donde £ es el módulo de Young. En función de la tensión elástica o las fuerzas a través de las dos caras de la rodaja son 
"= —G AU y F! =0'AÚg. 


Cuando se está propagando una perturbación, la tensión elástica (1) sobre cada sección no es la misma por lo que una 
rodaja de espesor dx estará sometida a una fuerza neta 
do OF 
F!'- F=dF = A—dz = —dz 
Ox Ox 
hacia la derecha, donde se ha tomado como positivo el sentido del sistema de fuerzas que da una tracción sobre la rodaja. 
Como la masa de la rodaja es dm = pAdz y su aceleración 9%u/0t? la ecuación de movimiento de la rodaja es 


do u 00 uu 
adds Ay = 5 =p. 1-55 
dF = adm > A rad pÁdx Ta _ dr =P ( ) 
A partir de la ley de Hooke (1-54) podemos escribir 
=p (1-56) 
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y derivando respecto de x se tiene 
do, 0%u 
dx Ox? 
que sustituyendo en la ecuación anterior (1-55) nos da 


u Eu ¿0u (1257) 
—7 == > =Ccó E 
02 poz? Ox? 

ecuación similar a la ecuación de onda, que nos indica que el campo de desplazamientos u se propaga a lo largo de la 
barra con una velocidad 


c= e (1-58) 
P 


El campo de tensiones elásticas gd satisface una ecuación similar, 


o Edo A [o 


927 pon a pda: 
así como el campo F, 

PF ESF ¿EF 

dE ii 


Propagación de la energía 


Si consideramos los distintos casos de propagación de ondas mecánicas vemos que en todos ellos se tiene un movimiento 
de partículas alrededor de una posición de equilibrio. Estas partículas del medio en el que se propaga la perturbación no 
tienen un movimiento promedio fuera de su posición de equilibrio. Como ya hemos dicho anteriormente, no se propaga 
la materia en sí, sino su estado de movimiento, lo que implica que lo que se está propagando con la onda es cantidad de 
movimiento y energía. 


Ya estudiamos anteriormente el caso de una cuerda tensa (1-19). Consideremos ahora el caso de las ondas elásticas longi- 
tudinales que se propagan a través de una barra. En un corte arbitrario de la barra, la fuerza que ejerce el lado izquierdo 
sobre el derecho es —0 A, hacia la izquierda si se tiene una tracción (con d positivo), produciendo un desplazamiento du 
y realizándose un trabajo —o Adu. Para que las ondas mecánicas comiencen a propagarse por la barra y que esta propaga- 
ción se mantenga es necesario comunicar de forma permanente a un extremo de la barra una potencia. Ésta es la potencia 
que se empleará en producir y difundir las ondas a lo largo de la barra. Para ello en nuestro ejemplo se comunicará una 
energía al lado izquierdo de la barra, que se propagará a lo largo de la misma. La potencia o trabajo por unidad de tiempo 
que el lado izquierdo transmite hacia el derecho será de la forma 
0W cos Qu 


Ot 0t 


(1-61) 


Cuando la perturbación pasa de una sección a otra, se transmite esa potencia. Para que se produzca propagación de 
izquierda a derecha debe suministrarse energía al extremo izquierdo de la barra. Esta perturbación puede ser aislada, con 
un pulso de extensión limitada producido por un suministro de energía en un instante dado, o continua, para lo cual se 
necesita un suministro continuo de energía. 


Energía que se propaga con una onda armónica 
Un caso particular muy interesante lo constituye una perturbación armónica, que da lugar a una onda armónica de la forma 
u = uysen(kx — wt). Sustituyendo en la ecuación (1-61) y derivando se tiene 


e = —0A [—wWuy cos(kx — wt)). 


Asimismo, de la ecuación (1-56) se tiene 


du : 
=E=—=E os(ka — 
g de kuo cos(kx — wt) 
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que junto a las relaciones w = ke y c = y E/p nos da 


o0W 

ide CA [pw*uz cos*(kz — wt)) (1-62) 
que es una expresión variable, pero siempre positiva. Observemos que esta expresión es de la forma genérica f(kzx — wt), 
por lo que satisface la ecuación de ondas. Esta onda corresponde a una onda de energía que se propaga con la perturbación, 
Normalmente no nos interesará la potencia instantánea que se transmite, sino la potencia media en un intervalo de tiempo, 
dada por 


(5 )- CA [pu < cos*(kz — wt) >). 


Ahora bien, el valor promedio de cos? en un periodo es 1/2, por lo que la expresión anterior queda como 


MA Sl 22 
(5 )> cAf zoo 18). (1-63) 


La energía en el medio en el que se propaga la onda mecánica será en parte energía potencial de deformación del medio 
(p.ej., compresión) y en parte energía cinética asociada al movimiento de las partículas del medio. En un cierto elemento 
de volumen, las proporciones de una y otra variarán, pero la energía total permanecerá constante si no hay disipación de 
energía. 

Observemos que la ecuación (2-7) de la referencia [4] nos da la energía total de un oscilador libre como E, = ¿mu A?, 
Particularizando a nuestro caso con una amplitud uy, una frecuencia angular w y tomando energías por unidad de volumen 
(con lo que en vez de la masa m tendremos la densidad p), se tiene 


z S 
2,2 
Ef = ze u0 (1-64) 
para la energía por unidad de volumen o densidad de energía en la barra a causa de las oscilaciones del movimiento 
ondulatorio. Si sustituimos esta ecuación en la (1-63) se tiene 
o0W 
—=r )=cA€;. (1-65) 
Ot 
Como c es la velocidad de propagación, c£;, es el flujo promedio de energía por unidad de área y de tiempo, mientras 
que al multiplicar por la sección A se tiene la energía promedio que fluye por unidad de tiempo a través de una sección 
transversal de la barra como consecuencia del movimiento ondulatorio. 


Se define la intensidad de la onda como el flujo promedio de energía por unidad de área y de tiempo 
1 W 
I==2 (Er). (1-66) 


Se obtienen resultados similares en un fluido para las ondas planas de presión y para las ondas esféricas lejos de la fuente, 
como se verá más adelante. 


1.7. Propagación de ondas mecánicas en fluidos. 


Las ondas de presión que se propagan en un fluido compresible reflejan variaciones de presión y densidad producidas por 
el movimiento de las moléculas en la dirección longitudinal en torno a su posición de equilibrio, como se verá en la figura 
1-6. En este desplazamiento no se produce un transporte neto de materia, aunque se transportará energía y cantidad de 
movimiento al igual que en una barra fina. 


Las diferencia de presiones entre distintos puntos actuará como fuerza restauradora sobre las moléculas del fluido. Como 
los fluidos no tienen resistencia elástica (fuerza recuperadora) a los cortantes, aunque sí tengan resistencia viscosa, en un 
fluido pueden propagarse únicamente ondas longitudinales en las que el movimiento del medio tiene lugar en la dirección 
de propagación. 
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1.7.1. La ecuación de ondas para un fluido no viscoso 


Cuando estudiamos la propagación de ondas de presión en un fluido hay que tener en cuenta que el fluido sin perturbar 
ya estará sometido a una cierta presión py en el equilibrio. Así, la presión total en un punto del fluido será la suma de la 
presión en el equilibrio y el cambio de presión g producido por la perturbación, 


p=po+p (1-67) 


La deducción de la ecuación de ondas para fluidos que vamos a presentar está basada en la bien conocida ecuación de 
continuidad para fluidos, 


> 0) 
Vp + =(pU) =0 (1-68) 
Ot 
y en la ley de Hooke para la compresión hidrostática, 
—p = KO (1-69) 


donde 9 es la deformación unitaria de volumen 9V/V y XK el módulo de compresibilidad. A menudo se escribe la ley de 
Hooke anterior para la compresión hidrostática (1-69) en función de la condensación 


9p SV 
IO (1-70) 
Po V 
en la forma 
p = Ko (171) 


Si consideramos la deformación de un elemento prismático infinitesimal de fluido con dimensiones dx, dy y dz a lo largo 
de los ejes coordenados, como consecuencia de la perturbación estas dimensiones pasarán a ser du + du, dy + duy y 
dz + du,, en función de las componentes del vector de desplazamientos al igual que vimos en el caso de las ondas en 
sólidos en la sección 1.6.1. Del mismo modo, la deformación unitaria de volumen queda 


dV OU 0Uy  0u% >, 
- > = Y 1-42 
“SE A e + dy le dz u ( ) 


0 


igual a la divergencia del vector desplazamiento en ese punto. Salvo en lo que respecta a la perturbación que se está 
propagando por el mismo, el fluido está en reposo, con po uniforme, por lo que Y; == Ve. Asimismo consideraremos 
que el término predominante en la densidad es la densidad en el equilibrio y que la perturbación no produce una excesiva 
variación en la densidad, por lo que podemos considerarla aproximadamente constante en la derivada temporal que aparece 
en la ecuación de continuidad (1-68). Así, ésta queda 

00 0% 
9 = VPF 


Aplicando el operador nabla a ambos lados de la ecuación de continuidad se obtiene 


Ve + po =0. (1-72) 


Sustituyendo 0 por el valor dado por la ley de Hooke (1-69) se tiene 


o? p 
Vo =p e, 
ESPERA 
que da la ecuación de ondas para un fluido, 


— = a = av? (1-73) 
0 


ds E (1-74) 
PO 


es la velocidad de propagación de las ondas en el fluido. 


donde 
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1.7.2. Energía de un fluido sometido a una onda de presión 


Estudiemos ahora el cambio en la energía de un elemento de un fluido por el que pasa una onda de presión. Dicho 
elemento, que se encuentra en 7 y tiene un volumen dV sufrirá en cada instante un desplazamiento u al paso de la onda. 
La energía cinética de este elemento vendrá dada por 


1 dy? 1 duy?* 


donde u es el módulo de u. Por otra parte, partiendo de la definición de módulo y coeficiente de compresibilidad para un 
cierto proceso, se tiene 
1 1 0Vv Vo 
A E > dV = -——dgp 
K Yo Op En 
para ese proceso!, Así, la energía potencial del elemento de volumen infinitesimal sería el trabajo realizado por p : 0 = p 
y queda como 


dV f? gp? p? 
dE, =-— dV = =— dp = Td = ¡5 0V 176 
"l y EL pap = y ao med ( ) 
ya que K = poc? de acuerdo con (1-74). La energía total instantánea por unidad de volumen del fluido será 
dEin 1 (duY.  g 
=% Es E . 1-77 
dVo ga ( di ) s 2poc? ( ) 


En esta expresión ú y g son funciones de 7 y t y por tanto la densidad de energía por unidad de volumen también lo será. 
Aunque esta expresión es válida para cualquier onda de presión en el fluido, para cada caso concreto debe conocerse la 
forma de ú y g. Más adelante veremos cómo es esta relación para los casos particulares de ondas planas y esféricas de 
presión sin disipación de energía. 


1.7.3. Velocidad del sonido en los fluidos 
Gases 


Examinemos las consecuencias de la ecuación (1-74) para la velocidad de las ondas de presión en una columna de gas. 
En esta sección veremos como se puede relacionar esta presión con otras propiedades de los gases. La ecuación de estado 
de un gas real, a bajas presiones o altas temperaturas, y en general lejos de las zonas en las que se produce licuefacción 
puede expresarse en una buena aproximación mediante la ecuación de estado de los gases ideales. En función de ésta, la 
densidad se puede escribir como 
Moi 5 m Mp 

pV =nRT = =RT > =pp,T)= 5 = 55 
en función de la masa molecular M, de la constante de los gases RR, de la presión p y de la temperatura 7”. La temperatura 
del gas tenderá a subir en aquellas zonas que se comprimen al paso de la onda y a bajar en aquellas que se expanden. 
Dependiendo del posible intercambio de calor entre las zonas calientes y frías se tienen dos posibles aproximaciones. 


Aproximación isoterma o de Newton. Esta aproximación fue propuesta por Newton en 1686 en base a consideraciones 
teóricas antes incluso de que se hiciera una distinción clara entre los procesos adiabáticos e isotermos. Si el proceso 
de propagación de la onda fuera muy lento, la temperatura en las zonas perturbadas por el movimiento ondulatorio se 
mantendría básicamente constante, al tener el sistema perturbado tiempo suficiente para alcanzar el equilibrio térmico con 
su entorno. En este caso se tendría un proceso isotermo, dado por la ecuación pV = cte y el módulo de compresibilidad 
sería 

p=cteV"! > K =p, 


propagándose la perturbación con una velocidad 


(1-78) 
P 


TK y x deben calcularse para el proceso adecuado. Aceptando la aproximación de Laplace del apartado 1.7.3 ésto sería para un proceso adiabático. 
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quedando la dependencia con la temperatura en la forma 


c=4/ E T 
M 


En la práctica esta aproximación da valores de la velocidad del sonido en el aire considerablemente menores que los 
valores experimentales. 


Aproximación adiabática o de Laplace. Fue Laplace quien en 1816 señaló que este proceso debe ser considerado 
adiabático. Aunque en realidad hay un pequeño flujo de calor, como éste da lugar a una atenuación de la onda y el hecho 
cierto es que el sonido puede viajar grandes distancias con una pequeña atenuación (aparte de la debida al factor 1/77), 
la magnitud del flujo térmico debe ser muy pequeña y por tanto la aproximación adiabática se presenta como una buena 
aproximación. 


Veamos las cosas de otra forma, el movimiento ondulatorio asociado a la propagación del sonido en los gases da lugar a 
perturbaciones locales que son lo suficientemente rápidas para que el intercambio de calor que tiene lugar pueda conside- 
rarse nulo, es decir, que el proceso es básicamente adiabático. Ésto será así mientras que la longitud de onda sea grande 
comparada con el recorrido libre medio de las partículas del aire. La absorción debiera aumentar para longitudes de onda 
más próximas al recorrido libre medio, pero en la práctica las longitudes de onda del sonido son mucho mayores que el 
recorrido libre medio (en factores grandes, de cerca de un millón). Asimismo, aunque el proceso es rápido y por tanto 
irreversible se considera en una primera aproximación como reversible con una ecuación de la adiabática pV? == cte, con 
y el índice adiabático, de forma que el módulo de compresibilidad X- vendría dado por 


H 


y Op 
> E Vo =D 


eós 5 (1-79) 


O RN (1-80) 


p=ctev"" 


y la velocidad del sonido en un gas sería 


que es mucho más próxima a los valores experimentales. En las tablas 1-1 y 1-2 se presentan distintos valores de la 
velocidad del sonido y de otros parámetros físicos en algunas sustancias. 


En realidad siempre hay un pequeño intercambio de calor, relacionado con una disipación de energía de la onda y por 
tanto con una atenuación de la misma. Como en la práctica el sonido puede viajar grandes distancias con una pequeña 
atenuación, el flujo de calor debe ser pequeño y por tanto la aproximación adiabática es esencialmente correcta. 


Salvo en el límite de bajas presiones y/o altas temperaturas los gases reales siguen ecuaciones de estado más complejas 
que la ecuación de estado del gas perfecto. Esto hace que el índice adiabático y sea en general función de presión y 
temperatura y que las expresiones anteriores sean más complejas. 


Líquidos 


La velocidad del sonido en los líquidos depende del módulo de compresibilidad y de la densidad al igual que en los gases. 
Sin embargo no es posible hacer un tratamiento simplificado como el que se ha hecho en los gases. En su lugar se utilizan 
diversas relaciones empíricas. 


1.7.4. Ondas planas longitudinales de presión en una columna de gas. 


Estudiaremos ahora las soluciones de la ecuación de onda que representan ondas planas de presión en un fluido. Ondas 
de este tipo pueden producirse moviendo un pistón en un extremo de un tubo, despreciando el efecto del rozamiento, 
del intercambio de calor y de la pequeña deformación de las paredes. Así, consideremos la columna de gas de la figura 
a lo largo de la cual se propaga una perturbación. Análogamente al caso de la barra que vimos en una sección anterior, 
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Figura 1-6 Propagación de una onda de presión en una columna de gas. Las distintas partículas de gas a lo largo de la columna oscilan alrededor de 
su posición de equilibrio. Esto da lugar a regiones en las que la presión aumenta (en tono más oscuro) y regiones en las que la presión disminuye (en 
tono más claro). La primera figura muestra la posición original de una serie de superficies equiespaciadas en la columna de gas antes de la perturbación, 
mientras que en las figuras inferiores se muestra cómo se ha desplazado cada una de las superficies y la forma de la onda de desplazamientos. la onda 
de presiones y la de velocidades están desfasadas 1/2 radianes respecto de ésta. Las líneas verticales de trazos muestran, sobre las otras figuras las 
posiciones originales de cada superficie mostradas en la figura superior. 


Figura 1-7 Propagación de una perturbación a lo largo de una columna de gas 
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Tabla 1-1 Velocidad del sonido y otras propiedades para algunos gases 


Gas c(ms71) p(kg/m79) poc(kgm"?s71) y 
Hidrógeno (H») 1270 0,090 114 1,41 
Nitrógeno (Na) 337 1,25 421 1,40 
Aire 331 1,29 428 1,40 
Oxígeno (O2) L/ 1,43 453 1,40 
Anhídrido carbónico (CO) 258 1,96 508 1,30 
Tabla 1-2 Velocidad del sonido y otras propiedades para algunos líquidos 

Líquido c(ms7!) p(kg/m"%)  poc10'(kg m”?s7!) 

Agua (4%C) 1418 1000 1,42 

Agua (25%C) 1493 997,1 1,49 

Mercurio 1450 13600 19,70 

Etanol 1210 790 0,96 

Metanol 1130 790 0,89 


estudiaremos cómo afecta la perturbación a una rodaja cilíndrica de espesor infinitesimal dx. Cuando no hay perturbación 
la rodaja está a una distancia x del origen de coordenadas y su espesor es dx. Cuando el sistema está perturbado, la rodaja 
se desplaza a 1 + u y su espesor pasa de ser dz a ser du + du. 


Como en este caso la propagación es unidimensional en la dirección x, el frente de onda es plano y la ecuación (1-73) 
pasa a ser 


e Pl: (1-81) 
con una velocidad de propagación c = y K/po. El desplazamiento satisface una relación similar, 


du Ku 
7 ES 2 077 (1-82) 


que en ambos casos corresponde a una onda plana que se propaga sin deformación en la dirección x con una velocidad c. 


El sonido no es más que una serie de ondas de presión que se propagan en el aire. En el desarrollo que hemos efectuado 
se ha considerado el caso más habitual en el que la presión varía relativamente poco respecto de su valor de equilibrio. Si 
esta variación fuera muy alta, como puede ser el caso de una explosión cerca del origen de la misma, las aproximaciones 
que hemos utilizado no serían válidas y el cálculo matemático se complica. De hecho ni siquiera la ley de Hooke lineal 
(1-69) sería válida en esas condiciones. 


Amplitud de desplazamiento, velocidad y presión de una onda armónica plana de desplazamiento. Impedancia de 
la onda 


Consideremos una onda plana armónica que se desplaza en la dirección x. El desplazamiento expresado en notación 
compleja será de la forma 


¡(ku—wt 
u= AS 0% [u = uy cos(kx — wt)] 
mientras que la velocidad sería 
e. 08 . (ha—wt 
v=ú= q —iwu E" E (u = uguw sen(kx — wt)] 


y la presión en forma compleja, sustituyendo en la ecuación (1-69) y teniendo en cuenta (1-42) 


du 


p==k0= ke Ko = —iku Ke 790 lo = Kkuo sen(ka — wt)]. (1-83) 


lo que muestra que la presión oscila en torno a su valor promedio con una amplitud po = Kkug. Como K = poe?, la 
expresión anterior queda 
2 
po = c“pokuo 
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o, eliminando k mediante la relación k =w/e 
fo = CPWUY = 2TCPOVUA (1-84) 


Observemos también que si el desplazamiento va como un coseno, la diferencia de presión oscila como un seno, esto es, 
cuando el desplazamiento alcanza sus valores extremos, la diferencia de presión es nula y viceversa. 


La impedancia de la onda (también llamada impedancia acústica específica) es el cociente entre la diferencia de presión 
y la velocidad. Para este caso en que ambas están en fase, queda como el cociente entre las amplitudes correspondientes, 


Es KK e kpoc? 
w w 


teniendo en cuenta que w = kc. El hecho de que la impedancia sea real indica que para una onda plana la presión acústica 


y la velocidad de las partículas del fluido que se mueven en la dirección de propagación están en fase. 


Propagación de la energía para una onda plana de presión en un fluido 


Para las ondas planas de presión en un fluido, la onda de desplazamiento es de la forma u = uy cos(kx — wt) y la onda 
de presión de la forma dada por (1-83), por lo que (1-77) da 


d En 1 21,2 
e dd ES + 23 | uy sen” (kx — wt) = pow”uy sen*(ka — wt), (1-86) 
0 0 


donde hemos utilizado ke = w y K = poc? y que recuerda a la expresión (1-64) para la energía de un oscilador libre. El 
promedio temporal de esa expresión da 


2 
Ez, e — Pes 
2poc? poc? 


1 
E = zPow*uj = 53 (cpowug) (1-87) 


2p0C 
La potencia promedio y la intensidad se obtienen de forma análoga a como se obtenían en (1-65) y (1-66) quedando, para 
la intensidad 


E Pes 174 Pes (1-88) 
2p0c poc 2% Z 


que es constante a lo largo del desplazamiento de la onda plana si no hay disipación de energía. En realidad siempre 
hay una disipación de energía en esta propagación, produciéndose calor, y la intensidad de la onda disminuye durante la 
propagación de la misma de forma aproximadamente exponencial, 


Ma-ES. aaa (1-89) 
donde xy es un punto de referencia, I(xp) la intensidad en dicho punto, a: la constante de atenuación del medio, de 
unidades? [neper/m]. En realidad, es ma = 2a: la magnitud que se utiliza habitualmente en acústica de salas y que 


representa la constante de atenuación del medio para la intensidad, también con unidades [neper/m]. La atenuación de la 
amplitud eficaz de la onda de presión tiene lugar de la forma 


p(a) = p(a0)€ 7 0=%) (1-90) 


en función de «+ que como vemos representa la constante de atenuación del medio para la amplitud de la presión eficaz. 


El neper es una unidad de medida del nivel de una magnitud o de atenuación, es similar al belio, pero definida en función del logaritmo neperiano, 


L(Np) = In (7) ; L(B) = log (5) 


donde Fy es el valor de referencia. 
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1.7.5. Ondas esféricas de presión en un fluido. 
La ecuación de ondas para una onda esférica que se propaga en un fluido y la solución para la onda de presión. 


Consideremos ahora las ecuaciones para una onda de simetría esférica, es decir, una onda que se propaga de forma radial, 
de forma que la presión es función únicamente de r y £, siendo r la distancia desde el origen. En este caso la ecuación de 
ondas (1-73), expresada en coordenadas esféricas, es de la forma 


%p 209 _10 (22) _ 10 


= 0 (1-91) 


Ñ Or 


Ti 


al tener únicamente dependencia radial. Si se realiza el cambio de variable g(r,t) = p(r, t)/r, la ecuación anterior queda 
en una forma más sencilla, 


o 10% 

Or c? 0t 
que tiene el mismo aspecto que la ecuación de ondas correspondiente a una onda plana que se propaga sin deformación 
en una dimensión con velocidad c, y que sabemos tiene la solución genérica 


p(r,t) = pi(r —ct) + palr + ct) 


compuesta de ondas que se propagan en una u otra dirección, obviamente con la misma velocidad que se obtuvo para las 
ondas planas en un fluido, 


as de (1-74) 
PO 
por lo que la onda real de presión será de la forma genérica 


1 1 
p(r,t) =p— po = ¿bl — 0) + rba(r + ct) (1-93) 


con una dependencia con la inversa de r junto a la de £. El primer término representa ondas esféricas salientes, que se 
alejan del origen y cuya amplitud disminuye de forma inversamente proporcional a la distancia al centro de la esfera, 
mientras que el segundo término representa ondas esféricas entrantes, que convergen hacia el origen y cuya amplitud 
aumenta con 1/r según se acercan. 


Presión, desplazamiento y velocidad local para una onda armónica esférica. 


Estudiemos ahora en mayor detalle una onda de presión armónica saliente. Ésta tiene la forma 
p=p-— po = golr) cos(kr — wt) = da cos(kr — wt) (1-94) 


donde la amplitud de la onda de presión es po = do/r (con do que no depende de 7) que disminuye con la distancia a la 
fuente. Esta misma ecuación se podría haber escrito en forma exponencial compleja, siendo (1-94) su parte real, 


pa ye e ae (1-95) 


Supondremos asimismo que la onda se propaga indefinidamente (o es completamente absorbida en las superficies a las 
que llega, como ocurre en una cámara anecoica). 


Para obtener la forma de la onda de desplazamientos y la velocidad local para una onda esférica partiremos de la ecuación 
de continuidad 


"E PRA 
Vp+ ¿e =0 (1-68) 


El operador gradiente tiene, en coordenadas esféricas, la forma 


y-(92109_1 3 
— XOr' r00” rsend0 0p 
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por lo que en simetría esférica Ye = 0g/0r queda como 


ó dp. 1 pi(krw ik pikr— Z 
Ve = Li, = do |-3€ E co! d. 
Or y r 


Si consideramos que la densidad varía poco respecto de la densidad en el equilibrio py se tendría, a partir de la ecuación 
de continuidad (1-68), 


0v 0% = 0 Po 1 i(kr—wt) 
=P = =P = = ly = — |-— + ik Uy. 
Ao SAA A | é is 
Integrando respecto del tiempo se obtendría la velocidad local, 
, 0 lt ¡(kr=wt) E Po |1 i(kr=wt) 
—=po00 = — |-+k e dr > u=-—|-+k|l Ur, (1-96) 
rw |r powr |r 


e integrando de nuevo, el desplazamiento de las partículas del fluido al paso de la onda, 


1 ¡(kr—wt) -, mn 1 z 7 —wt) -, 
pi Y ES a [Eat] a (1-97) 
rw? |r powér | r 


A grandes distancias de la fuente el término en 1/7 se hace despreciable y el desplazamiento puede expresarse con buena 
aproximación como 


¿= —D et 


pow?2r 
cuya parte real es 
k 
eos a sen(kr — wt) = go(r) sen(kr — wt) = uy(r) sen(kr — wt) (1-98) 
powér powc 


Como se ve, para puntos alejados del emisor, la amplitud de la onda de desplazamientos también disminuye como 1/r, 
de forma inversamente proporcional a la distancia a la fuente, 


u =up(r) sen(kr — wt) = E sen(kr — wt). 
donde la relación entre uy y (o viene dada, a grandes distancias, por 
polr) = cpopwun[r) = 2rcpyvuo[r) + do = cpow(o = 2rrcpovCo, (1-99) 
expresión similar a la (1-84) que obtuvimos para las ondas planas, pero con uy y fo dependiendo de r. Hay que insistir 


en que ahora estas expresiones son válidas únicamente a grandes distancias de la fuente. 


Impedancia de una onda esférica en un fluido. 


La impedancia de la onda esférica en un punto del fluido queda como 


q=*- rF Po/r 5 PE =- E $1 (1-100) 
¿/kr 
[544] 00/0007 + 
r r 
que consta de una parte real y de una parte imaginaria. Lejos de la fuente la parte imaginaria se hace cero y la parte real 


se puede aproximar por 
Z = poc (1-101) 


lo que significa que lejos de la fuente de nuevo ambas magnitudes, presión diferencial y velocidad, están en fase. 


Se suele denominar resistencia acústica específica a la parte real de la impedancia y reactancia acústica específica a la 
parte imaginaria. 
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Energía de una onda esférica lejos de la fuente. Intensidad de la misma. 


Para obtener la intensidad de la onda esférica (1-94) lejos de la fuente utilizaremos la recién obtenida ecuación (1-98) y 
la ecuación (1-77) teniendo en cuenta que ahora la amplitud es uy(r) = Co/r en vez de uy. Así, la energía por unidad de 
volumen o densidad de energía queda como 

Es loa 2(q) = 1 pow*c5(r) _golr) Pd 


A geo “o 2 y2 —2c2p e 20? por? 


que disminuye como 1/r?. Puede verse que la relación entre la intensidad de una onda esférica y la densidad de energía 
es de la forma (1-66), por lo que se tiene para la intensidad lejos de la fuente emisora, 


pa(r) _ Peli) _ par) _ Perl) _ 0 do 


T=«c€; = Es e 
“— 2pye poc 2Z Z 2cpor2  r? 


(1-102) 


donde ] 
IS UI 7 
2cp Cp 2Z  Z 
la intensidad de una onda esférica es inversamente proporcional al cuadrado de la distancia a la fuente emisora. La energía 
que atraviesa en la unidad de tiempo una superficie esférica de radio r es por tanto, 


0 > (1-103) 


OW 
—- |= I 4rr? = 4r1o (1-104) 
Ot 

que ya no depende de r como era de esperar, ya que la energía total que atraviesa las distintas superficies esféricas debe 
ser la misma si no hay disipación de energía. Como en realidad siempre hay disipación de energía en forma de calor, la 
intensidad para una onda esférica disminuye de la forma aproximada 


—2a(r—ro 2 —2a(r—ro 2 —Malr—ry 
Ane 1(r) =1(r0) (2) esto - ra) (3) gormar=y (1-105) 


donde rg es un punto de referencia, 7 (ro) la intensidad en dicho punto, a: la constante de atenuación del medio para la 
presión [neper/m]. y mo la constante de atenuación del medio para la intensidad [neper/m]. Para la amplitud eficaz de la 
onda de presión se tiene, 


2 —QUU TO 
p(r) = pro) (2) Ext” (1-106) 


Energía en el caso general 


En el caso más general la potencia por unidad de superficie transmitida por la fuente al medio puede escribirse como 
W =úp =ú2Z (1-107) 


y consta de dos partes, una real y una imaginaria. La parte real avanza con la onda mientras que la parte imaginaria se 
emplea en dilatar la superficie de la onda, quedando acumulada en forma de energía potencial. Esta componente, a la que 
se denomina reactiva tiene un valor apreciable sólo en la primera fase de la propagación. 


1.7.6. Ondas superficiales en un líquido. 


Las ondas superficiales que se forman en la superficie libre de un líquido tienen una naturaleza distinta de las anteriores, 
ya que no son ni como las ondas de presión ni como las ondas en una membrana elástica tensa (la superficie libre del 
agua), sino bastante más complicadas. 


Al perturbar la superficie se generarán ondas, con alternativamente valles y crestas, pero ésto no puede corresponder a 
movimiento simple arriba y abajo, ya que el agua no puede ir para abajo al ser prácticamente incompresible, sino que 
debe desplazarse a los lados, saliendo de donde está. Lo que ocurre en la realidad es que ya no tenemos ni una onda 
transversal ni una longitudinal, sino una mezcla de ambas, ya que el agua que se encuentra cerca de la superficie se mueve 
aproximadamente en circunferencias, para el caso de ondas de longitud de onda larga en aguas profundas. Este estado de 
movimiento se aprecia en la superficie del agua como una ondulación que se propaga con una cierta velocidad. 


La velocidad de propagación estas ondas es mayor para aguas profundas que para aguas poco profundas, y ésta es la razón 
de que las olas rompan cuando llegan a la playa. El libro de Feynmamn[3] tiene una interesante discusión sobre ésto. 
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1.8. Ondas electromagnéticas. 


Hasta ahora hemos visto diversos ejemplos de ondas mecánicas longitudinales y transversales que se propagan a lo largo 
de un medio. No necesariamente el movimiento ondulatorio implica un movimiento mecánico, sino que comprende cual- 
quier tipo de perturbación que se propaga a través de un medio. Por ejemplo una perturbación de un campo eléctrico o 
magnético que se propaga a través de un medio también entra dentro del movimiento ondulatorio. Este es el caso de las 
ondas electromagnéticas, en las que un campo eléctrico y un campo magnético armónicos, perpendiculares entre sí y a la 
dirección de propagación se propagan en el espacio, como se ve en la figura 1-8. 


Figura 1-8 Una onda electromagnética está formada por una perturbación armónica del campo eléctrico y del campo magnético que se propagan en el 
espacio con velocidad c. Ambos campos están en fase y sus amplitudes están relacionadas por E = CB. 


1.8.1. Ondas planas 


Consideremos las ecuaciones de Maxwell del electromagnetismo para el vacío (no hay cargas libres ni corrientes, p = 0 
y j= = ()), para una solución particular consistente en un campo eléctrico E según el eje y y un campo magnético B según 
el eje z perpendiculares entre sí, y a la dirección de propagación x, 


E,=0E,=EE,=0 => E=(0,E,0) 
B¿=0B,=0B,=B >= B=(0,0,B) 


quedando las ecuaciones de Maxwell como 


E 
> =0 Ley de Gauss para el campo eléctrico (1-108a) 
OB 2 
e > 0 Ley de Gauss para el campo magnético (1-108b) 
OE e 
E 0 Ley de Faraday-Henry (2) (1-108c) 
Se =- de Ley de Faraday-Henry (x) (1-108d) 
5 =() Ley de Ampéere-Maxwell (y) (1-108€e) 
OB 
DER Si = Eno S Ley de Ampére-Maxwell (1) (1-108f) 


Como se deduce de las ecuaciones (1-108a), (1-108b), (1-108c) y (1-108e), ni E ni B dependen de y o z. Por tanto los 
campos E y B dependen sólo de x y £, teniendo en cada instante cualquiera de ellos el mismo valor en cualquier plano 
perpendicular al eje x. 
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Si derivamos la ecuación (1-108d) respecto de x se tiene 


E a 0B 
0x2 0x0t 

y derivamos (1-108f) respecto de t se tiene 
0?B 0%E 


Oxrot = €oMo ve 
que combinadas dan 
e _ 1 0E 
02 — opa Ox? 
que es de la forma de la ecuación de ondas (1-2), lo que indica que la solución anterior corresponde a un campo africo 
que se propaga con una velocidad 3 


(1-109) 


1 


= —== (1110) 
v/coHo 
como una onda plana, y que se puede expresar en la forma 
E = E(1-—ct) (1111) 
Para el campo magnético también se obtiene la ecuación de ondas en la forma 
9? il PB 
OB_ 10B (1112) 
0t2  eouo Ox? 
lo que implica que el campo magnético también es de la forma 
B= B(x-ct) (1113) 


Igualmente se podrían haber utilizado las soluciones E(x + ct) y B(x-+ct), pero por simplicidad utilizaremos únicamente 
las z — ct. 


1.8.2, Ondas planas armónicas 


Consideremos ahora el caso particular de ondas planas armónicas de frecuencia y = w/27 y longitud de onda A = 27 /k. 
En este caso particular la forma de los campos eléctrico y magnético es, para propagación hacia la parte positiva del eje 
de las x, 

E = Ep sen k(x — ct) = Ep sen(kz — wt) 

B = Bosen k(z — ct) = Bosen(kz — wt) 


donde hemos utilizado la relación w = kc. Estas expresiones deben cumplir la ecuación (1-108d), teniéndose para cada 
derivada parcial 


(1-114) 


OE 
— = kEpcos(kx —wt 
da kEo cos(kx — wt) 
OB 
> —wBy cos(kx — wt) 


que al sustituir da la relación entre las amplitudes del campo eléctrico y magnético, 
Eo = E Bo =cBo >  E=cB (1-15) 


que implica directamente que los campos E y B están en fase, ya que c es una magnitud real y por tanto Ly y Bo lo están 
(no hay ningún factor complejo que introduzca un desfase) y lo mismo E y B que tienen factores similares. 


La que acabamos de ver no es la única solución aceptable de las ecuaciones de Maxwell en forma de onda plana que se 
propaga según la dirección x (que satisface además la ecuación de ondas para propagación según 2), ya que hay otras que 
satisfacen la condición de que el vector resultante del producto vectorial E x B esté dirigido en el sentido de propagación 
de la onda electromagnética, que es la que da las relaciones de signos correctas para que se satisfagan las ecuaciones de 
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Maxwell. Por ejemplo, y manteniendo la dirección de propagación a lo largo del eje x con sentido positivo, una solución 
con el campo eléctrico dirigido según +2 y el campo magnético según —y también sería válida. 


Otra solución posible es aquélla en la que en todo momento los campos eléctrico y magnético se mantienen perpendicula- 
res entre sí y a la dirección de propagación, que permanece fija, pero rotando el plano de oscilación de los campos eléctrico 
y magnético a lo largo de la propagación en torno a ésta. En este caso se tiene una onda polarizada circularmente, mientras 
que en el caso anterior en el que permanece invariable el plano de oscilación de los campos eléctrico y magnético se tiene 
una onda polarizada linealmente. En el caso de una onda circularmente polarizada, los campos eléctrico y magnético se 
expresarían de la forma 

E, = Eo sen(kx —wt) E, = +E0cos(kzx — wt) de 

B, = FBocos(kzx — wt) B, = Bosen(kz — wt) : 
donde hay dos posibilidades de signos dependiendo de si el giro es a derechas o a izquierdas. 


1.8.3. Energía de una onda electromagnética 


En los textos de electromagnetismo (Véase por ejemplo [2]) se ve que la densidad de energía asociada con el campo 
eléctrico de una onda electromagnética en el vacío es 


1 
En = ¿co? (1-117) 
y que la densidad de energía magnética en el vacío es 
1 
En = rd (1-118) 
que teniendo en cuenta que ego = 1/c? y que E = cB queda 
1 1 1 
AA E? = -60E?. 
2 20 2poc? qe 


la densidad de energía eléctrica de una onda electromagnética es igual a la densidad de energía magnética. Por tanto la 
densidad de energía total del campo electromagnético es 


Eorm = eoE* (1119) 


1.8.4. Intensidad de una onda electromagnética 


La intensidad instantánea de la onda electromagnética (energía que atraviesa la unidad de área en la unidad de tiempo) 
será 
I(x,t) = €ogme= ceo E? (1-120) 
mientras que la intensidad promedio será 
=—— 1 
I=Eogmc= ¿eco Es (1-121) 


al tener E forma armónica. 


1.8.5. El vector de Poynting 


Consideremos ahora el producto vectorial c?ey E x B para el caso de una onda electromagnética plana. Este vector apunta 
hacia la dirección y sentido de propagación y tiene como módulo 


I2eE x B| = e EB = ce E? 


que es la intensidad instantánea dada por la ecuación (1-120). Al vector 


de ExB (1-122) 
se le denomina vector de Poynting. El flujo de este vector de Poynting a través de una superficie A, 
= =p dE 
il “e BExBdA= (1-123) 
A dt 


es la energía que fluye a través del área A por unidad de tiempo, dada por otra parte por dE/d£. 
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1.9, Velocidad de grupo y velocidad de fase 


Aunque ya hemos hablado de ondas que pueden descomponerse en serie de ondas armónicas no nos hemos detenido en 
las propiedades de la onda resultante, sino únicamente en las de cada una de las ondas armónicas en que se descompone. 
Cada una de estas ondas se propaga con una velocidad dada por la ecuación (1-8) e = w/k, que representa la velocidad 
a la que se propaga una cierta fase, o un cierto estado de la perturbación, y denominamos a esta velocidad velocidad de 
fase de la onda. La velocidad de propagación de la onda no necesariamente coincide con la velocidad de movimiento de 
ninguna entidad física, por ejemplo la velocidad de propagación de una onda plana de presión en el aire no coincide con 
la velocidad de las moléculas del aire por las que pasa la onda, solamente representa la velocidad a la que se propaga el 
estado de movimiento de las mismas, no su velocidad. 


En principio una onda armónica “pura” no podría transmitir más información que su frecuencia y su longitud de onda, pero 
en la práctica no es posible hablar de una onda armónica “pura” más que para analizar su propagación en un desarrollo de 
Fourier, ya que cualquier onda debe tener un principio y un final cuando llega a un receptor. Esto hace que en esta onda 
con principio y fin tengamos en realidad una serie de ondas superpuestas que a través de una serie o integral de Fourier dan 
la onda que estamos estudiando. De hecho aquí solo podríamos hablar de velocidad de fase en el intervalo de existencia 
de la onda. Por tanto, para transmitir información real mediante una onda ésta no puede ser una onda armónica, sino que 
necesitaremos ondas no armónicas y en general no periódicas. 


En el caso más general, la velocidad de propagación de la onda dependerá de la frecuencia (o la longitud de onda) de la 
misma. En este caso decimos que se produce dispersión y se denomina medio dispersivo al medio en el que tiene lugar 
este comportamiento y medio no dispersivo en caso contrario. 


Como consecuencia de ésto, si por ejemplo las ondas sonoras de alta frecuencia viajaran más rápido que las de baja 
frecuencia, un pulso de sonido se oiría a una cierta distancia como una sucesión de sus distintas componentes, primero 
las de alta frecuencia y después las de baja frecuencia. Análogamente ocurriría con un pulso de luz blanca, que al viajar a 
través de un medio dispersivo se vería, si las bajas frecuencias viajaran más rápido que las bajas, primero como luz roja, 
luego verde y luego azul. La experiencia diaria nos dice que ésto no es así en el aire, en el que las ondas sonoras viajan de 
forma esencialmente independiente de la frecuencia y lo mismo ocurre con las ondas electromagnéticas viajando a través 
del aire o del vacío, pero puede ocurrir cuando estas ondas se propagan a través de distintos medios. Así, en general, 
únicamente las ondas estrictamente armónicas se propagarían conservando la forma. 


Cuando se tiene un medio dispersivo, si c disminuye con la frecuencia en algún rango de frecuencias, se dice que el medio 
presenta dispersión normal en ese rango mientras que si aumenta con la frecuencia en ese rango se dice que el medio 
presenta dispersión anómala en ese rango. Esta es una terminología que tiene su origen en la óptica. 


normal anómala a 


A E y 


Figura 1-9 Dispersión normal y anómala 


Consideremos ahora el caso de una onda de forma arbitraria que se propaga en un medio dispersivo. En medios lineales 
esta onda puede descomponerse en una serie ondas armónicas mediante el análisis de Fourier. La velocidad de fase de cada 
una de las componentes armónicas viene dada por la ecuación de onda para esa onda concreta. En un instante posterior, 
se puede obtener la forma de la onda resultante sumando las ondas individuales que la componen. Recordemos que si la 
onda es periódica, las componentes se obtienen como una serie de Fourier, mientras que si la onda es no periódica, las 
componentes se obtienen a través de una integral de Fourier. Ahora bien, si el medio es dispersivo, cada componente se 
propaga con distinta velocidad y por tanto, una onda de perfil arbitrario que se propaga en un medio dispersivo cambia de 
forma según se propaga. Denominaremos velocidad de grupo a la velocidad con que se propaga el pulso. 


Estudiaremos lo que queremos decir con ésto a partir del paquete más sencillo, el formado por dos únicas ondas armónicas 
de igual amplitud, pero frecuencias ligeramente distintas, con dw = wj — wa, por lo que sus longitudes de onda también 
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diferirán levemente, con 9k = k1 — ka 


Ape tera) (o+4w)] 


X= (1-124a) 


iNk—i e" 
Ls ¿en 50k)u—(w 30)] (1-124b) 


Así, estas ondas tienen la misma amplitud, pero los números de onda difieren en dk y las frecuencias angulares difieren 
en óww, con k el número de onda promedio y w la frecuencia angular promedio de las dos ondas. 
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Figura 1-10 Paquete de ondas resultante de la composición de dos ondas de la misma amplitud y frecuencias y números de onda próximos. El tiempo 
corre hacia abajo y la línea en la parte superior de cada gráfico representa el movimiento de la envolvente con la velocidad de grupo vy, mientras que la 
línea en la parte inferior representa el movimiento con la velocidad de fase del paquete Uf. 


Aunque en ciertos puntos en el medio las ondas estén en fase en un instante determinado, en otras posiciones o en otros 
instantes de tiempo esto no será así. En este caso, la onda resultante tiene un aspecto similar a las pulsaciones que se ven 
en la sección 1.4.1 y la figura 1-7, ambas de la referencia [4], pero con una particularidad nueva muy importante; como se 
muestra en la figura 1-10 tanto la envolvente como la onda resultante misma se están propagando, no están quietas en el 
espacio sino que se están propagando en el mismo y con velocidades en general distintas. Es precisamente a la velocidad 
con que se propaga la envolvente a lo que denominamos velocidad de grupo o velocidad con que se propaga el pulso. 


Aunque cada una de las ondas viaje con su velocidad de fase (ligeramente distinta para cada una de ellas si el medio es 
dispersivo e igual para ambas si no lo es) el patrón de pulsaciones viajará con la velocidad de grupo, que puede llegar 
a ser bastante distinta del promedio de las velocidades de fase de las ondas individuales, que será w/k, al ser k y w los 


respectivos valores promedio. Denominamos velocidad de fase del paquete a este promedio dado por el cociente 
w 
Cf=>, 1-125 
da ( ) 
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donde ahora w es la frecuencia angular promedio y k el vector de onda promedio, en vez de las frecuencias y vectores de 
onda de cada una de las ondas individuales. 


Analicemos un poco más la onda resultante. La superposición de las dos ondas puede escribirse como 
aaa iii EN ete Pl ei ko=ut) 
y la onda física por la parte real de esta ecuación, 
x =2Acos E o 74) cos(kx — wt) (1-126) 


La onda resultante tiene un número de onda k que es el promedio de los números de onda de las dos componentes y una 
longitud de onda A = 27 /k. Como se ve en la figura 1-10, la amplitud de la onda está modulada, teniendo la envolvente 
un número de onda 0k/2, la mitad de la diferencia entre los números de onda de las dos componentes. Aunque la longitud 
de onda de la modulación es 47 /8k, la separación entre dos pulsaciones consecutivas es la mitad de esa distancia, 277 /0k. 


Si en vez de considerar la evolución espacial de la elongación x de la onda en un determinado instante de tiempo, se 
considera la evolución temporal de la resultante en un punto determinado se tiene un comportamiento similar donde ahora 
en vez del eje x se tiene el eje de tiempos £. Las pulsaciones viajan entonces con la velocidad 


o) = d0 101 o =D 

2 "Sk  ki—ko “dk 

que es la velocidad a la que hemos denominado velocidad de grupo. Nótese que hemos utilizado el hecho de que las 
frecuencias de ls ondas son muy próximas. A partir de las relaciones w = ke y A= 27/k, ésta se escribe como 


(1-127) 


c=c+Hk=e- A (1-128) 


Si se conoce la dependencia c(k) de la velocidad de fase con el número de onda es posible calcular la velocidad de 
grupo del paquete. Para un medio no dispersivo c es independiente de k y por tanto cf = Cy, no habiendo en este caso 
diferencia entre las velocidades de grupo y de fase del paquete. Sin embargo, en un medio dispersivo sí habrá diferencia y 
la velocidad de grupo puede ser tanto menor como mayor que la velocidad de fase, dependiendo de si se tiene dispersión 
normal o anómala, y se producirá deformación del paquete. 


Hay otra cosa que debemos destacar, la velocidad de grupo representa la velocidad de transmisión de la señal únicamente 
en condiciones de dispersión normal o en general, cuando la velocidad de grupo es menor que la velocidad de fase. De 
hecho, jugando con valores arbitrarios de w y A en (1-127) se podría obtener una velocidad de grupo arbitrariamente 
grande a partir de ondas relativamente lentas, o incluso que la envolvente se propague en sentido contrario, lo que obvia- 
mente no puede corresponder a la velocidad de propagación de una señal. Otro caso es el de una onda electromagnética 
que se propaga en un medio dispersivo en condiciones de dispersión anómala, donde la velocidad de grupo puede ser 
mayor que la velocidad de la luz en el vacío, como ocurre con la propagación de rayos X en un cristal de carbono, Ésto es 
menos sorprendente de lo que parece, lo que no puede superar la velocidad de la luz en el vacío es la velocidad de trans- 
misión de una señal o la velocidad de movimiento de una partícula, y ésto no ocurre en ninguno de los casos anteriores. 
En realidad no es posible transmitir información a una velocidad mayor que la menor de las velocidades de fase de las 
ondas individuales. 


Aunque la ecuación (1-127) que se acaba de obtener es válida de forma estricta solo para el paquete tratado, que contiene 
dos ondas ligeramente de frecuencias y números de onda ligeramente distintos, sigue siendo útil para paquetes de ondas 
armónicas más complejos en los que interviene un rango de frecuencias, mientras que en ese rango w sea una función que 
varía lentamente con k. 


1.10. El efecto Doppler! 


Cuando la fuente y el observador se encuentran en movimiento relativo, la frecuencia con que llegan al observador es 
distinta de la frecuencia con que han sido emitidas por la fuente, debido al fenómeno conocido como efecto Doppler. En 
la figura 1-11 se muestra la causa de este efecto. Cuando la fuente emisora está en movimiento, los máximos de la onda 
emitida llegan con mayor frecuencia cuando la fuente va al encuentro del receptor, y con menor frecuencia cuando la 
fuente se aleja del receptor. 


tLlamado así en honor al científico austríaco Christian Doppler, 1803-1853 
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Figura 1-11 Efecto Doppler para una velocidad de movimiento de la fuente emisora menor que la velocidad de propagación de la onda. 


1.10.1. El efecto Doppler en una dimensión. 


FE — veót —| F : al 


Figura 1-12 Efecto Doppler en una dimensión 


Consideremos la situación que se muestra en la figura 1-12. Un emisor de ondas F se mueve hacia la derecha con una 
velocidad v¿ mientras que el receptor O se mueve también hacia la derecha con una velocidad v,. En el instante inicial 
t = 0 en el que el emisor emite una onda, la distancia entre el emisor y el receptor es lo, llegando la onda al receptor en 
un tiempo t. Durante ese tiempo el receptor ha recorrido una distancia v,t y la distancia total que la onda ha recorrido en 
ese tiempo t es ly + v,t, que si ces la velocidad de propagación de la onda se puede escribir como 


lo 


C— Ur 


ct = ly + upt => b= 


Transcurrido un tiempo dt el emisor ha avanzado una distancia v¿Ót y emite una nueva onda que llegará al receptor 
transcurrido un tiempo t' medido desde el origen de tiempos común. La distancia recorrida por la onda será ahora (ly — 
VeÓt) + vyt' que, ya que esta segunda onda ha viajado durante un tiempo (+t' — $t) a una velocidad c nos da 


, lo +(c— ve)ót 


c(t' — $t) = lo — veÓt + Yyt! > t 
C— Ur 


por lo que el intervalo en el que llegan al receptor las dos ondas emitidas por el emisor con una separación ót es 


Si el emisor emite ondas con una frecuencia y, en el intervalo de tiempo dt habrá emitido vót ondas. Como esas mismas 
ondas las recibe el receptor en un intervalo de tiempo dt”, la frecuencia con que el receptor reciba las ondas es y! = v0t/0t' 


(ya que vdt = 1'ót”), que se escribe como 
y Ct, 


(1-129) 


Pp = 
C—Ve 


y que nos da la relación entre la frecuencia y con que emite la fuente y la frecuencia v! con que el observador recibe la 
señal. 


Cuando la velocidad de desplazamiento de la fuente emisora v¿ y la velocidad de desplazamiento del receptor v, son 
pequeñas comparadas con la velocidad c de propagación de la onda, es posible simplificar la ecuación anterior aprove- 


chando el desarrollo en serie 
(Eye) 
C 0 
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quedando 


A OT O 
rei dead ir el ara 


donde se ha despreciado el término de segundo orden v,ve/c?. En función de la velocidad relativa del receptor respecto 
de la fuente emisora, V7. = U, — Ve, esta expresión se escribe como 


Ys (1 BA me) y (1-130a) 
Cc 
que, teniendo en cuenta que w = 21 y se puede expresar para la frecuencia angular, 
es (1 da 22) 0 (1-130b) 
Cc 


Cuando el observador se aleja de la fuente, v,. es positiva y la frecuencia medida es menor (la longitud de onda es mayor), 
mientras que si el observador se acerca a la fuente, 0,2 es negativa y la frecuencia medida es mayor (longitud de onda 
menor). 


1.10,2. El efecto Doppler en el espacio 


Cuando la fuente y el observador no se mueven en la misma dirección, sino de forma que sus direcciones de movimiento 
forman un ángulo 0, las ecuaciones anteriores quedan en la forma 


es (1 0 cos 0) y (1-131a) 
e 

ei = (1 e cos 0) W (1-131b) 
Cc 


1.10.3, Ondas de choque 


Figura 1-13 Cuando la velocidad de la fuente emisora es mayor que la velocidad de propagación de la onda se tienen ondas de choque. Se denomina 
cono de Mach a la envolvente de esas ondas de choque. 


Un caso especial se tiene cuando la velocidad relativa entre el observador y la fuente es mayor que la velocidad de 
propagación c de la onda. En este caso, la fuente avanza con más rapidez que el frente de onda. La superficie tangente a 
todas las ondas sucesivas es un cono cuyo eje es la línea recta descrita por la fuente en su movimiento y cuyo ángulo a de 
apertura viene dado por 

sena == (1-132) 

Ve 

En este caso el movimiento ondulatorio es una onda cónica que se propaga en las direcciones perpendiculares a la envol- 
vente, a la que se denomina onda de Mach u onda de choque. Esta onda es el sonido brusco que se escucha cuando 
un avión que viaja a una velocidad superior a la del sonido pasa cerca de nosotros, aunque también se observa en otras 
circunstancias, como cuando un barco de mueve sobre el agua a una velocidad mayor que la velocidad de las ondas 
superficiales en el agua. 
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Figura 1-14 Una bala atraviesa un fluido a una velocidad menor (izquierda), ligeramente mayor (centro, mach 1,01, la forma redondeada favorece una 
alta presión delante de la bala que hace que el cono esté ligeramente por delante) y claramente mayor (derecha) que la del sonido en el mismo. 


EY 


Figura 1-15 Un avión F/A-18 rompe la barrera del sonido. (Foto de la US Navy realizada por John Gay). 


1.11. Reflexión, refracción y difracción de una onda 


1.11.1. El principio de Huygens. 


Christiaan Huygens (1629-1695), astrónomo, matemático y físico holandés nacido en La Haya, propuso un método para 
estudiar la evolución del frente de onda de un movimiento ondulatorio. En su trabajo utilizó modelos geométricos elabo- 
rados sobre papel que comparaba con los resultados experimentales. 


Para ello partió de la suposición de que cada uno de los puntos del medio, perturbado al ser alcanzado por el frente de 
onda, actúa como un nuevo foco de perturbación que emite ondas esféricas con la misma velocidad, frecuencia y longitud 
de onda que la onda inicial. Todos los puntos del frente de onda emiten simultáneamente dichas ondas (denominadas 
ondas secundarias o elementales). 


El principio de Huygens establece que la envolvente de las ondas secundarias emitidas desde el frente de onda constituye 
el nuevo frente de onda. En todos los puntos que no están en un determinado instante sobre esa envolvente se produce una 
interferencia destructiva entre las distintas ondas secundarias que hace que su efecto se cancele. De esta forma podemos 
construir un frente de ondas conociendo el mismo frente en un instante anterior. 


1.11.2. Reflexión y refracción en una superficie plana de separación entre dos medios. 


Cuando las ondas que se propagan en un medio homogéneo llegan a una superficie donde las propiedades del medio 
cambian de forma abrupta, se producen los fenómenos de reflexión y refracción (o transmisión) de la onda. La onda 


34 Capítulo 1. Movimiento ondulatorio 


reflejada es la que se propaga desde la superficie de separación por el mismo medio (D) que la onda incidente mientras que 
la onda refractada es la que se propaga por el segundo medio. 


refractada 


incidente 


Figura 1-16 Reflexión y refracción en la superficie plana de separación entre dos medios O y Y. 


Consideremos la siguiente situación. En el fluido (D) una onda plana de frecuencia angular w y número de onda k¡ = w/c1 
incide sobre la superficie formando un ángulo 91 con la normal a la misma. Consideremos al plano de la figura, que es el 
plano de incidencia de la onda incidente. La onda reflejada, si existe, tendrá la misma frecuencia y por tanto un vector de 
onda de la misma magnitud que la onda incidente. Denominamos ángulo de incidencia al comprendido entre la dirección 
de propagación de la onda incidente y la normal a la superficie de separación entre dos medios en el punto de incidencia. 
Denominamos ángulo de reflexión al comprendido entre la dirección de propagación de la onda reflejada y la normal, y 
ángulo de refracción al formado por la dirección de la onda refractada y la normal. 


Se puede ver que, para la denominada reflexión regular, el ángulo de incidencia debe ser igual al de reflexión, y ambos 
rayos deben estar en el mismo plano, 


'=01. (1-133) 


Hemos hablado de reflexión regular ya que en la reflexión hay que tener en cuenta la relación entre la longitud de onda 
A de la onda incidente y el tamaño de las irregularidades de la superficie. Si la longitud de onda de la onda incidente es 
menor que el tamaño de las imperfecciones se tiene reflexión regular y la onda se refleja con un ángulo de reflexión igual 
al de incidencia y poca distorsión. Si por el contrario la longitud de onda de la onda incidente es del orden o mayor que 
el tamaño de las irregularidades de la superficie se tiene la denominada reflexión difusa, en la que se produce reflexión en 
un rango más amplio de ángulos y distorsión de la onda reflejada respecto de la incidente. 


Las leyes de la reflexión son muy importantes para el estudio de la acústica de locales, y a menudo se colocan las 
superficies de forma que las reflexiones intensifiquen el sonido en ciertas zonas. Obviamente, la menor dimensión de estas 
superficies, a las que se denomina tornavoces, debe ser mayor que la longitud de onda del sonido incidente para que no se 
produzca distorsión del mismo. 


También puede verse que para la refracción debe cumplirse la relación 
k¡ sen 0, = ko sen 0, 
también con ambos rayos en el mismo plano. Esta relación se escribe habitualmente en la forma 


sen 0; ko C1 
ss E, 1-134 
sen 02 ki C2 ( ) 


Esta relación se conoce como ley de Snell, especialmente cuando se trabaja con luz. En ese caso se dice que el coeficiente 
c1/c2 es el índice de refracción del segundo medio respecto del primero. En Óptica se suele utilizar el término índice 
de refracción de un medio de forma absoluta, refiriéndolo siempre a la velocidad de propagación de la luz en el vacío, 
n = C/Cmedio. Dependiendo de la velocidad de propagación de la onda en ambos medios, el ángulo de refracción será 
más o menos próximo a la normal que el ángulo de incidencia. 


Como los medios por los que se propaga una onda no son en general completamente homogéneos, la refracción es un 
fenómeno muy frecuente. A causa de esta inhomogeneidad las propiedades del medio cambian de un punto a otro y por 
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tanto la velocidad de propagación de las ondas en ellos. Por ejemplo, la velocidad de propagación del sonido en el aire 
cambia de un punto a otro debido a que su temperatura no es uniforme y por tanto se produce refracción. Así, en un día 
soleado las capas de aire próximas a la superficie terrestre se encuentran más calientes que las más altas y por tanto la 
velocidad del sonido, que aumenta con la temperatura, es mayor en las capas bajas que en las altas. Ello da lugar a que 
el sonido, como consecuencia de la refracción, se desvíe hacia arriba. En estas condiciones, la comunicación entre dos 
personas suficientemente separadas se ve dificultada. Por las noches tiene lugar el fenómeno contrario, ya que la Tierra se 
enfría más rápidamente que el aire. 


La refracción de la luz es la causa de que el Sol parezca más grande al amanecer que más tarde. El aire que está cerca de la 
superficie terrestre es más denso que el estrato de aire superior, y así sucesivamente, hasta que la capa de aire desaparece 
a unos 400 kilómetros. En estas condiciones, la luz emitida por el borde superior del Sol atraviesa la atmósfera con un 
ángulo distinto que la luz del borde inferior. Esto da lugar a dos ángulos de difracción divergentes, que producen la ilusión 
óptica que nos hace creer que el Sol es más grande. La refracción es también la causa de que los cuerpos sumergidos en 
agua tengan un tamaño aparente distinto del real. 


1.11.3. Difracción. 


El fenómeno de la difracción se presenta siempre que una onda se encuentra con un obstáculo, o una abertura, de dimensio- 
nes comparables o menores que su longitud de onda. El obstáculo puede ser una abertura, un muro o en general cualquier 
elemento que bloquee una parte del frente de onda (Véase la figura 1-17). 


Figura 1-17 Difracción de ondas en el agua a través de una rendija. 


En estas condiciones, la propagación del movimiento ondulatorio es rectilínea únicamente en la dirección no afectada por 
el obstáculo, pero aparece una propagación no rectilínea desde las direcciones correspondientes al borde del obstáculo, 
con una intensidad que depende de la dirección considerada. Hay que destacar que la magnitud de este efecto depende de 
la relación entre las dimensiones del obstáculo y la longitud de onda de la onda incidente. Si por ejemplo, estrechamos la 
abertura, la difracción se hace más marcada, los máximos laterales se separan del máximo central en forma de abanico, 
mientras que si ensanchamos la abertura, los máximos laterales se cierran sobre el centro y el efecto de la difracción va 
haciéndose progresivamente menor y la propagación a través de la abertura cada vez más rectilínea. Cuando el diámetro 
de la abertura es muy grande en comparación con la longitud de onda prácticamente no hay difracción y la propagación 
es esencialmente rectilínea 


Podemos entender la difracción a partir del principio de Huygens. Así, los bordes del obstáculo actuarían como origen 
de ondas secundarias. Estas ondas secundarias se propagan en todas las direcciones y dan la impresión de que la onda 
incidente ha rodeado el obstáculo. 


La difracción de las ondas sonoras permite que en determinadas ocasiones un sonido se oiga aunque el receptor no esté 
en línea recta respecto del emisor (y no haya reflexiones), a causa de obstáculos en el camino de la onda sonora. Para un 
obstáculo dado cuanto mayor sea la longitud de onda del sonido mayor será el efecto de la difracción y cuanto menor sea 
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ésta menor será el efecto de la difracción. Esto hace que, por ejemplo, las ondas sonoras de longitud de onda alta (sonidos 
agudos) sean más propensas a producir sombra acústica que las de baja longitud de onda (sonidos graves). 


Alta frecuencia 


5 


Baja frecuencia 


Figura 1-18 Difracción para altas y bajas frecuencias. Sombra acústica. 


Es importante tener en cuenta estos efectos a a la hora de diseñar barreras acústicas, ya que a bajas frecuencias, tendrá 
lugar más difracción en los bordes de la barrera, con una menor sombra acústica y por tanto, con un peor apantallamiento 
por parte de la barrera. Por el contrario a altas frecuencias la difracción será menor y el apantallamiento mejor. 


Las ondas electromagnéticas visibles también experimentan los fenómenos de difracción. Sin embargo, como sus longitu- 
des de onda son muy pequeñas (400 nm a 700 nm, para el espectro visible), los efectos de difracción sólo serán apreciables 
en obstáculos y aberturas de dimensiones muy pequeñas. 


Sin embargo, el efecto de la difracción es apreciable en otras zonas del espectro. Por ejemplo, la ondas de radio en las 
frecuencias correspondientes a la denominada onda media (AM) tienen una longitud de onda grande comparada con el 
tamaño de la mayor parte de los obstáculos que encuentran en su camino, por lo que se produce difracción con cierta 
facilidad alrededor de muchos edificios y alcanzan lugares que no alcanzan ondas de radio de longitud de onda más corta. 


La ondas de radio en frecuencia modulada (FM) tienen una longitud de onda menor y por tanto la difracción con los 
edificios es menor, por lo que se reciben con menor calidad. Por ello en muchos lugares se recibe mejor la onda media 
que la frecuencia modulada. Este efecto es aún más marcado para las ondas de televisión cuya longitud de onda es mucho 
menor, por lo que es necesario poner las antenas en los tejados. En este caso el efecto de la difracción es positivo para la 
recepción. 


El efecto de la difracción también limita la resolución máxima de la instrumentación que utiliza reflexión para observar 
objetos pequeños, ya que se producirá reflexión regular únicamente si la longitud de onda de la onda utilizada es pequeña 
comparada con el tamaño del objeto. En caso contrario se producirá reflexión difusa (causada por la difracción) que impide 
resolver detalles de tamaño menor que la longitud de onda de la onda utilizada. Por ejemplo, no es posible, utilizando 
sonido audible, resolver algo con una precisión mayor de +2 cm ni observar átomos con luz visible. 


Da 


Capítulo 2 
Ondas estacionarias 


2,1. Ondas estacionarias en una cuerda 


2.1.1. Condiciones en los extremos de una cuerda 


Consideremos ahora la siguiente situación, un tren de ondas avanza a lo largo de una cuerda tensa y llega al extremo de la 
misma. Cada sacudida que llega al soporte ejerce una fuerza sobre éste, de forma que la reacción a esta fuerza actúa sobre 
la cuerda y da lugar a una sacudida reflejada que se propaga en sentido contrario. Sin embargo, esta sacudida reflejada 
debe siempre satisfacer una condición evidente, si el extremo está fijo la sacudida reflejada debe ser tal que éste debe 
permanecer en reposo. 


También puede darse el caso de que el extremo esté libre. En este caso, la sacudida que llega continúa propagándose, 
pero en sentido contrario, dando lugar también a una onda reflejada. Se denominan condiciones en los extremos a las 
condiciones que deben satisfacerse en los extremos de la cuerda. 


Cuando se estudia la propagación de esta perturbación en el tiempo, y en particular cómo se produce la reflexión en el 
extremo se ve que, cuando el extremo está fijo la sacudida se refleja, invirtiéndose los valores de la elongación y de la 
velocidad. Por el contrario, cuando la reflexión tiene lugar en un extremo libre se invierte la dirección de la velocidad, 
pero no cambia la de la elongación. 


Figura 2-1 En la figura se muestra la perturbación real y las perturbaciones virtuales con que se compondrá para satisfacer las condiciones en el extremo. 


Se puede considerar este proceso como si ocurriera de la siguiente forma. Supongamos que la cuerda se prolongara 
indefinidamente más allá del extremo, de forma que la perturbación continuara en la parte imaginaria de la cuerda, como 
si no hubiera condición alguna en el extremo. Al mismo tiempo, consideremos que una perturbación virtual, que ha estado 
propagándose en la parte imaginaria, se mueve hacia la parte real y forma la perturbación reflejada. Las características 
de esta perturbación virtual dependen de si se tiene un extremo fijo o libre, mostrándose los dos casos en la figura 2-1. 
Siempre podemos imaginar estas reflexiones suponiendo que la onda real y la onda virtual que sale desde el otro lado del 
extremo. 


Así, se tiene en el eje negativo de las = una perturbación real que se desplaza hacia el origen de coordenadas, del tipo 
u(x — ct) como corresponde a una onda que se desplaza hacia la derecha. Consideremos en primer lugar la perturbación 
virtual inferior, correspondiente al caso de extremo fijo. Esta perturbación es la imagen de la perturbación original pero 
reflejada dos veces, una primera respecto del plano yz y una segunda respecto del plano yx. Matemáticamente, ésto 
implicaría una inversión de la dirección x, es decir  —> —< y un cambio de signo de la perturbación, quedando ésta como 
—u(—zx — ct). Esta función es del tipo u(x + ct) (ya que es —u(—(x + ct))) correspondiendo por tanto a una perturbación 
que se desplaza hacia la izquierda con velocidad c. La perturbación virtual superior, correspondiente al caso de extremo 
libre, es simplemente la imagen de la perturbación original reflejada respecto del plano yz. Matemáticamente se obtendría 
igual que la perturbación virtual anterior, pero sin cambiar de signo la perturbación, únicamente invirtiendo el signo de 
la coordenada x. Así, sería una perturbación del tipo u(—x — ct) que corresponde a una onda que se desplaza hacia la 
izquierda con velocidad c (igual que antes, es u(—(% + ct))). 
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Figura 2-2 Reflexión de una perturbación en un extremo fijo (izquierda) y en un extremo libre (derecha). Las perturbaciones real y virtual se muestran 
con línea de trazos mientras que la perturbación resultante (superposición de las dos anteriores) se muestra con trazo continuo. La zona virtual está 
marcada en gris. 


En un instante cualquiera, la perturbación total en un punto se obtiene como la suma de los valores en ese punto de las 
perturbaciones original y reflejada. Este hecho se conoce como principio de superposición y es en general cierto siempre 
que las perturbaciones sean suficientemente pequeñas, y en el caso de elementos elásticos, que no se sobrepasen los límites 
elásticos. En la figura 2-2 se muestran las perturbaciones incidente, reflejada y total tanto para el caso de reflexión en un 
extremo fijo (parte izquierda de la figura) como para reflexión en un extremo libre (parte derecha de la figura). Como se 
ve en la figura, en el caso de reflexión en un extremo fijo las perturbaciones real y virtual se componen de forma que su 
suma, la perturbación resultante, es siempre nula en el extremo fijo. 


2.1.2. Ondas estacionarias en una cuerda fija a un extremo 


En vez de la propagación y reflexión de una perturbación aislada, consideremos ahora que ocurre cuando se tiene un tren 
de ondas armónicas que incide de forma continua sobre el extremo fijo de una cuerda. En este caso parece que en dicho 
extremo se produjera un tren continuo de ondas propagándose en el sentido contrario y siendo, en el rango de validez del 
principio de superposición, la perturbación resultante en cada punto la suma de las perturbaciones incidente y reflejada. 
En general denominaremos interferencia al efecto producido por distintos trenes de ondas que atraviesan simultáneamente 
una cierta región del espacio. 


En este caso se tiene un comportamiento peculiar, ya que el aspecto de la perturbación resultante no muestra ondas 
recorriendo la cuerda en sentidos opuestos. En vez de eso, si la frecuencia es suficientemente grande, la cuerda aparece 
subdividida en un cierto número de segmentos. La forma de la cuerda sigue siendo la de una función armónica pero, 
mientras que en una onda armónica que se propaga la amplitud de la onda permanece constante según avanza la onda, 
aquí la onda resultante permanece fija en posición, mientras que la amplitud varía como consecuencia de la oscilación. 
Existen unos puntos denominados nodos que permanecen siempre en reposo, mientras que los puntos equidistantes entre 
dos nodos consecutivos tienen amplitud máxima (denominándose a estos puntos vientres o antinodos). Se denomina a la 
vibración resultante en este caso onda estacionaria u onda confinada, y su característica fundamental es que se trata de 
una onda que oscila, pero sin una propagación neta, aunque sea el resultado de la superposición de dos ondas viajeras. 
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Figura 2-3 Onda estacionaria resultante de la superposición de dos ondas que se reflejan en un extremo fijo. 


Estudiemos analíticamente estas ondas estacionarias. Para ello consideremos el siguiente tren de ondas, 


0 


2 
u¡ = Y cos(kz — wt) = ul cos le — ct) 


que representa un tren de ondas armónicas incidente que se propaga hacia la derecha con velocidad c. Si el extremo 
derecho está fijo, la ecuación de la onda reflejada vendrá dada, de acuerdo con lo visto en la sección 2.1.1, por 


2 
ur = —ul cos(—kz — wt) = —ul cos o — ct) 
que puede escribirse como 
2 
uy = —u/ cos(kz + wt) = —ul cos la + ct) 


La onda resultante vendrá dada por la superposición de u; y Ur, 


2 
u =u; + uy = Y [cos(kx — wt) — cos(kx + wt)] = cos Ta — ct) — cos y + ct) 


Desarrollando los cosenos de la suma y diferencia de dos ángulos se obtiene 
0 0 217 
u = 2u sen kz sen wt = 2u” sen qe sen 2rrvt (2-1) 


En un instante de tiempo genérico (en el que sen wt no sea nulo), la forma (2-1) de la onda resultante en el espacio es una 
sinusoide [240 sen w+] sen kx cuya amplitud 2u" sen wt varía a su vez de forma sinusoidal con el tiempo. Para instantes 
de tiempo en los que se satisface la relación 


wi=2rvt =0,71,27,... 
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la amplitud es nula y la cuerda permanece recta. En estos instantes de tiempo las ondas incidente y reflejada están en 
oposición de fase y ambas ondas se cancelan. Por el contrario, para instantes de tiempo en los que se cumpla la relación 


wit =2rvt = 1/2,371/2,57/2,... 


la amplitud tiene un máximo de valor 2u%, que es la suma de las amplitudes de las ondas que interfieren. En esos instantes 
de tiempo, las ondas incidente y reflejada están en fase y dan por tanto, la máxima amplitud resultante. En los demás 
instantes, la amplitud tiene valores intermedios que varían de forma sinusoidal con el tiempo con una frecuencia angular 
w. 


¿Cómo oscila cada punto de la cuerda y cómo varía de un punto a otro de la cuerda el desplazamiento causado por la onda 
en un instante de tiempo dado? A partir de la ecuación (2-1) podemos representar esta oscilación de cada punto como un 
movimiento armónico de amplitud 2u/ sen kz y frecuencia . 


Así, a partir de la ecuación (2-1) (igual que a partir de la figura 2-3), se ve claramente que cada punto de la cuerda efectúa 
un movimiento armónico simple de la misma frecuencia y para todos los puntos, pero con una amplitud 2u" sen kz que 
depende de la posición del punto. Al contrario de lo que ocurre con una onda que se propaga, todos los puntos situados 
entre dos nodos oscilan en fase. 


Para valores de x que cumplen la relación 


271 
ka === =0,1,27,+*+ =n7 = (2n)1/2 (Mi 3/0,:L, 2, o) 
la amplitud de la oscilación es nula. Estos puntos permanecen siempre en reposo y son por tanto nodos. Por el contrario 
en los puntos donde se satisface la relación 


—kx = == =71/2,37/2,57/2,--- = (2n + 15 


(=0,, Zo.) 
la amplitud alcanza su valor máximo 2u0, siendo por tanto vientres o antinodos. Los nodos se encontrarán en aquellos 
puntos en los que se cumpla 


A 
zx =0,4/2,24/2,:-: e (2n)A/4 (n= 0,1; 25%) (2-2) 
y los vientres en aquéllos en los que se cumpla 
=x = A/4,31/4,51/4,---= 2n+1Í (n =0,1,2,..) (2-3) 


Los nodos consecutivos están separados entre sí en media longitud de onda, lo mismo que los vientres consecutivos. 


Todo esto se muestra de forma gráfica en la figura 2-3. Entre dos nodos sucesivos cualesquiera cada punto oscila sinusoi- 
dalmente con la misma frecuencia, pero el movimiento resultante no se propaga a otros puntos, sino que permanece fijo 
en el espacio. Ésta es la característica de los modos de vibración. Si nos encontramos con un movimiento en el que cada 
punto del móvil se mueve sinusoidalmente en torno a su posición de equilibrio con la misma frecuencia (aunque cada 
punto lo haga con distinta amplitud) tenemos un modo de vibración. 


2.1.3. Ondas estacionarias en una cuerda fija a ambos extremos 


Consideremos ahora el caso más habitual, una cuerda fija que se encuentra sujeta a ambos extremos. Un tren de ondas 
armónicas se refleja en dichos extremos y al estar ambos fijos deben ser los dos nodos, con un desplazamiento nulo de la 
cuerda en cada uno de esos extremos. Como los nodos consecutivos están separados entre sí por media longitud de onda, 
si la longitud de onda de las ondas armónicas es A, para que se cumpla la condición de desplazamiento nulo en los nodos, 
la longitud de la cuerda debe ser un número entero de semilongitudes de onda, 


ALALA A 
L==2,27,37,"=nNn=z A 2 
32393 de” (n 1,2 ) (2-4) 
Dicho de otra forma, si consideramos una cuerda de una determinada longitud L, sólo pueden producirse en la cuerda 
ondas estacionarias de determinadas frecuencias, aquéllas cuyas longitudes de onda son 21, 2L/2,2L/3,..., es decir, 


AS —, con (n= 1,2,..) (2-5) 


2.2. Ondas estacionarias en una columna de aire 41 


y distancias entre nodos L, L/2, L/3,.... Como la velocidad de propagación es la misma para todas las frecuencias, las 
permitidas serán 
e € 6 E 
= —>712>7 33357," =n>> = LL) — 
dde e PAN 


Al igual que en la descomposición de Fourier, se denomina a la frecuencia más baja vo = c/2L frecuencia fundamental 
o armónico principal mientras que a las frecuencias superiores se les denomina armónicos superiores (segundo armónico 
219, tercer armónico 34, etc.). 


Estas frecuencias son características del sistema oscilante y de la naturaleza de los límites. En general cualquiera de los 
movimientos posibles de la cuerda puede analizarse considerando la superposición de ondas estacionarias de distintas 
frecuencias, o dicho de otra forma, de una serie de modos de vibración, por supuesto siempre que estemos considerando 
un sistema lineal en el que es aplicable el principio de superposición. 


Hay una diferencia importante a resaltar entre un sistema de un muelle con un peso y una cuerda vibrante. Mientras que 
el sistema formado por el muelle y el peso tiene una única frecuencia angular natural de vibración dada por w¿ = k,/m, 
la cuerda vibrante tiene un número infinito de frecuencias propias, la fundamental y todos los armónicos. Si desplazamos 
de su posición de equilibrio un peso sujeto por un resorte este sistema oscilará a una única frecuencia de vibración. Si 
deformamos inicialmente una cuerda de modo que su forma coincida con uno de los posibles armónicos, cuando se libere 
la cuerda, ésta vibrará con la frecuencia correspondiente a dicho armónico en particular. 


¿Qué ocurrirá cuando la forma inicial de la cuerda no coincide con la de ningún armónico? Éste es el caso que se tiene 
cuando se golpea o pulsa una cuerda de forma arbitraria. En este caso, la perturbación inicial puede expresarse mediante un 
desarrollo de Fourier en función de ondas armónicas de longitudes de onda coincidentes con las de los posibles armónicos 
de la cuerda. De esta forma, la vibración resultante es una mezcla de ondas estacionarias, cada una de ellas correspondiente 
a una componente de Fourier de la perturbación inicial. Por ejemplo, al golpear una cuerda de un piano, estarán presentes 
la frecuencia fundamental y los distintos sobretonos. 


Por la misma razón, hay una diferencia adicional entre el muelle y la cuerda vibrante en relación a la resonancia. Mientras 
que el sistema peso-muelle puede entrar en resonancia únicamente a su frecuencia fundamental, una cuerda vibrante puede 
entrar en resonancia en cualquiera de sus armónicos. Para un sólido hay aún más posibilidades, ya que además de ondas 
estacionarias transversales como las que hemos visto en el caso de una cuerda, también puede haber ondas estacionarias 
transversales de torsión como las que aparecieron en el puente de Tacoma Narrows que se muestra en la figura 2-4, así 
como ondas estacionarias longitudinales. 


Las amplitudes de los distintos armónicos dependen de la forma particular en que se haya pulsado la cuerda inicialmente. 
El tono del sonido principal emitido por la cuerda está determinado principalmente por la frecuencia fundamental mientras 
que el timbre depende del número de armónicos que se hayan excitado y de sus amplitudes. 


Si se tuviera uno de los extremos de la cuerda libre, éste sería un vientre en vez de un nodo. 


Como hemos visto, la frecuencia fundamental de la cuerda vibrante tensa viene dada por 1) = c/2L, donde según se vio 
en (1-14) c = Y F/p, con F la fuerza tensora y p, la densidad lineal de masa. Así, la frecuencia fundamental de una 
cuerda vibrante de densidad lineal de masa p, tensa con una fuerza F es 


a (27) 
2L Y pi 

En los instrumentos de cuerda se encuentran excelentes ejemplos de aplicación de esta ecuación. Así, éstos se afinan 
variando la fuerza tensora F, de forma que un aumento de F aumenta la frecuencia y viceversa. Asimismo, cuanto mayor 
es la longitud de la cuerda, menor es la frecuencia. Por ejemplo, en un piano, las cuerdas correspondientes a los tonos altos 
tienen menor longitud que las cuerdas correspondientes a los tonos bajos, lo mismo que las cuerdas del violonchelo son 
más largas que las cuerdas de un violín y por tanto dan un tono más grave. Para obtener una baja frecuencia sin necesidad 
de tener F muy bajas o longitudes muy largas, tanto en los pianos como en las guitarras se recurre a forrar las cuerdas 
que dan las notas bajas con un hilo metálico, o a utilizar cuerdas más gruesas. 


2.2. Ondas estacionarias en una columna de aire 


Hemos estudiado en temas anteriores la propagación de una onda de presión a lo largo del gas contenido en un tubo. Si éste 
tiene longitud finita, las ondas se reflejarán en los extremos del tubo de la misma forma que lo hacen las ondas transversales 
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Figura 2-4 Un viento persistente excitó en el puente de Tacoma Narrows una onda estacionaria de torsión como la que se ve en la figura superior. Esta 
onda estacionaria entró en resonancia con el viento, provocando finalmente la rotura del puente el 7 de noviembre de 1940, apenas 4 meses después de 
su inaguración. La amplitud de la onda estacionaria generada fue de unos 4, 3 m, la frecuencia de 14 ciclos por minuto y su longitud de 853 m, siendo la 
densidad lineal de masa del puente de 6700 kg/m, llegando a acumular una energía de 4 x 107 J, Ésta se liberó rápidamente durante la rotura llegando 
a ser la potencia de 9 x 10% W durante los segundos en los que se produjo el colapso del puente. (Imágenes de dominio público de la Stillman Fires 
Collection: Tacoma Fire Dept.) 


que se propagan en una cuerda tensa cuando llegan a sus extremos. Así, tendremos también ondas estacionarias y el tubo 
con aire tendrá también un número infinito de frecuencias propias. 


Los tubos sonoros pueden ser de dos tipos, abiertos y cerrados. Son abiertos si tienen dos o más orificios y cerrados si 
tienen un único orificio. La mayoría de los instrumentos musicales de aire están construidos con tubos abiertos. Solo unos 
pocos instrumentos como algunos tubos de órgano utilizan tubos cerrados. 


Si el extremo en el que se produce la reflexión es cerrado, el desplazamiento longitudinal de las partículas en dicho 
extremo debe ser nulo y por tanto el extremo cerrado corresponde a un nodo. Ahora bien, ésto corresponde a un punto de 
variación máxima de presión, es decir a un vientre de presión, como se aprecia en la figura 2-5. 


Cuando el extremo del tubo está abierto las cosas son más complicadas y dependerán de la relación entre el diámetro 
del tubo y la longitud de onda de la perturbación. Cuando el tubo es estrecho comparado con la longitud de onda de la 
onda que se propaga en el mismo, como ocurre con la mayor parte de los instrumentos musicales, el extremo abierto se 
mantiene a la presión atmosférica y por tanto es un nodo de presión y un vientre de elongación (Un tratamiento riguroso 
muestra que, en realidad, el nodo está ligeramente fuera del tubo, a unas seis décimas de su radio). 


De esta forma, las ondas longitudinales de presión en tubo se reflejan en los extremos del mismo, abierto o cerrado, del 
mismo modo que las ondas estacionarias en una cuerda tensa lo hacían en los extremos fijo o libre. 


El efecto de los agujeros intermedios en un instrumento de viento es actuar como un extremo abierto y dividir la columna 
de aire en segmentos, produciéndose en cada uno de ellos una frecuencia propia. Así, tapando unos u otros agujeros, se 
tiene un nodo de presión (y un vientre de elongación) en cada agujero destapado, de forma parecida a lo que ocurre con 
el nodo de desplazamiento que aparece al pisar con el dedo la cuerda de un instrumento de cuerda. En general, en un 
instrumento de aire la longitud efectiva de la columna está menos definida que la longitud de una cuerda tensa fija por sus 
extremos. 


La excitación de la columna de aire se produce mediante una embocadura, de distintas formas según los instrumentos. Sin 
una embocadura adecuada, una corriente de aire proyectada a través de un estrecho canal tendería a romperse a su salida 
formando turbulencias. Sin embargo, si dicha corriente de aire se hace chocar contra una arista más o menos afilada (bisel 
en al caso de la flauta) o se hace vibrar una lengieta como en el caso de un clarinete, en vez de turbulencias se tienen 
oscilaciones periódicas y el sistema corriente-vibrador se convierte en un excitador. Así, la misión de esta embocadura es 
transmitir la vibración de la misma al aire que se encuentra en el tubo y de todas las frecuencias que ésta comunica las 
únicas frecuencias que sobreviven en la columna de aire son aquéllas que coinciden con la del armónico fundamental y 
los armónicos superiores, produciéndose resonancia en éstas. Para que se efectúe transmisión de la vibración la abertura 
en la que se encuentra la embocadura no puede ser un nodo de desplazamiento, pero no hace falta que sea un vientre sino 
que el punto de excitación puede estar en un lugar intermedio. Tampoco es necesario que las aberturas del tubo coincidan 
con los extremos del tubo, que pueden estar cerrados y tener el tubo varias aberturas a lo largo del mismo 
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Figura 2-5 Diagramas de varias ondas estacionarias en el aire contenido en un tubo abierto o cerrado. En la figura se representa la envolvente de la onda 
de presión. Cada uno de los nodos de la presión correspondería a un vientre de desplazamiento y cada uno de sus vientres a un nodo de desplazamiento. 


En la figura 2-5 se muestran la frecuencia fundamental y los dos primeros armónicos en tubos abiertos y cerrados, señalán- 
dose la envolvente de las ondas de presión en cada caso. Como se ve, para una misma longitud, la frecuencia fundamental 
de un tubo abierto es doble que la de un tubo cerrado. También se observa que las frecuencias de los armónicos para un 
tubo abierto van como 2, 3, 4 veces la frecuencia fundamental, mientras que los armónicos de un tubo cerrado van como 
3, 5,7 veces la fundamental. Así, en las vibraciones en un tubo abierto están presentes todos los armónicos mientras que 
en un tubo cerrado sólo lo están los impares y por tanto el timbre de ambos será distinto, aunque se hayan construido para 
que tengan la misma frecuencia fundamental. 


La transmisión de las vibraciones a la columna de gas también puede hacerse a través de un extremo móvil como en el 
tubo de Kundt. Éste consiste en un tubo provisto de una válvula para la entrada de gas, cerrado en sus extremos por dos 
émbolos, uno de ellos ajustado a un dispositivo vibrador (que puede ser simplemente el frotamiento de una varilla ajustada 
al mismo) y otro de ellos de posición ajustable. Cuando se espolvorea en el interior del tubo un polvo fino y se excita 
la vibración del primer émbolo, ésta se transmite al gas y ajustando convenientemente la distancia entre los émbolos se 
seleccionan las ondas estacionarias producidas, que se pueden apreciar mediante los montones que forma el polvo. Cada 
montón representa un nodo de desplazamiento, ya que el polvo ha sido desplazado de todos los demás lugares en los que 
el gas se movía al nodo, donde no lo hace. La distancia entre dos nodos sucesivos dará la mitad de la longitud de onda del 
modo de vibración. 


2.3. Ondas estacionarias en una membrana tensa 


Analizaremos ahora las ondas estacionarias que se forman en una superficie, con sus correspondientes modos normales. 
Para ello utilizaremos los resultados que obtuvimos en la sección 1.5 para una membrana tensa y particularizaremos a dos 
casos, el de una membrana rectangular y el de una membrana circular, ambas fijas en sus bordes. 
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2.3.1. Ondas estacionarias en una membrana rectangular fija en sus bordes 


Como sugiere la simetría del problema, partiremos de la solución (2-16) e impondremos las condiciones de contorno en 
los límites de la membrana, u = 0 para [x = 0, a] y [y = 0, b] si la membrana tiene dimensiones a Xx b. Una solución del 
tipo 

u(x,y,t) = Asen k, x sen kyeteo (2-8) 


se hace 0 para x = 0 y y = 0, mientras que si se cumple además que tanto k.¿a como kb son un número entero de veces 
Tr la condición de desplazamiento nulo también se satisface en los otros dos bordes. Así, si ka = lr y kyb = mr, con 
l,m = 1,2,3,... las frecuencias de los distintos modos normales de vibración vendrán dadas por 


2 


w l my? 
k= La 4/62 +2 => w(l,m) = re (5) + (5) . (Q-9) 


13 


Figura 2-6 Modos de vibración de una membrana rectangular con una dimensión el doble de la otra (b = 2a). Para cada modo se muestra el par (1, m) 
correspondiente. Los signos + y — muestran la fase relativa del movimiento de las distintas partes de la membrana. 


La frecuencia fundamental correspondería a m = l = 1 y los sobretonos con m = l serán armónicos de la frecuencia 
fundamental, pero en general aquéllos con m X l no lo serán. En el caso más general, ni unas frecuencias serán múltiplos 
de las otras ni serán múltiplos de ningún número salvo el caso especial mencionado. Si a? y b? son inconmensurables 
(es decir, a?/b? no puede expresarse como el cociente de dos enteros), todas las frecuencias serán distintas. En caso 
contrario pueden existir dos o más frecuencias idénticas para diferentes pares de índices m, l; se dice en este caso que se 
trata de modos de vibración degenerados. Para un modo de vibración, las líneas nodales en la membrana constan de dos 
conjuntos de rectas dirigidos cada uno de ellos según uno de los ejes coordenados, y dividen la membrana rectangular en 
Lx m rectángulos. 


La vibración más general de una membrana rectangular puede representarse como una superposición lineal de los distintos 
modos de vibración, en la forma 


[0,0] 00 
l ; SS 
u(z, y,t) = Ey Ey Aim Sen =— sen - grat) (2-10) 


donde Wim Y Am deben determinarse a partir de las características geométricas de la membrana y de las condiciones 
iniciales de cada problema concreto. 


2.3.2. Ondas estacionarias en una membrana circular de radio R fija en sus bordes 


Como vimos en la sección 1.5.3 si la membrana circular existe en el centro, no pueden aparecer las funciones de Bessel 
de segunda especie en la solución de la ecuación de ondas (1-23). Si el radio de la membrana es RR, los modos normales 
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Figura 2-7 Modos de vibración en una membrana circular fija en su borde. Para cada modo se muestra el par (m, n) correspondiente. Los signos + y 
— muestran la fase relativa del movimiento de las distintas partes de la membrana. 


de vibración deben corresponder a los valores de k que anulen la función de Bessel de primera especie Jn (k.R). Para 
cada uno de los ceros de la función de Bessel ,,, a los que denominaremos ¿nn refiriéndonos al cero enésimo, existirá 
el correspondiente k y por tanto, los números de onda de las ondas posibles en la membrana deben cumplir 


KEmnR = Ímn Q-11) 


y los modos normales de vibración han de ser de la forma 


Umnl[r,0,t) = AmnJm(kmnr) cos(mó — y EPA (2-12) 
con las frecuencias angulares naturales dadas por 
on = Zin = Kimi = PS. (2-13) 


El ángulo de fase azimutal ymn, que es una de las constantes que aparecen en esta solución, debe ser ajustado a partir de 
las condiciones iniciales, ya que determina las direcciones que pasan por el centro de la membrana a lo largo de las cuales 
se tendrán líneas nodales de desplazamiento nulo y ésto depende claramente de las condiciones de la excitación inicial. 


Por otra parte, si observamos cualquiera de las funciones de Bessel que aparecen en la figura 1-4, para el modo asociado 
al primer cero jmo no habrá líneas circulares internas en las que los puntos de la membrana no oscilen, al no haber 
ceros dentro. Sin embargo, cuando analizamos los modos correspondientes a los siguientes ceros vemos que además 
de en el borde de la membrana tampoco habrá oscilación en los puntos de la membrana correspondientes a uno de los 
ceros interiores. Como estamos en el mismo modo en el cero exterior y en los interiores, debe cumplirse kmn.Rk = Jmn 
para el cero exterior (que anula el desplazamiento en el borde de la membrana y caracteriza al modo) y kmnTi = JÍmi 
para cada uno de los ceros interiores que aparecen para cada nodo a una distancia 7, del centro de la membrana (con 


1=1,...,n-— 1). Como se trata del mismo modo, km» es el mismo en ambos casos y las distancias r;, que vienen dadas 
por 
n=bÉR «uu isleon=L Q-14) 
Imn 


definen una serie de puntos a lo largo de una circunferencia en los que la elongación es nula para el modo (m, n), a los 
que se denomina círculos nodales. 


Así, para un modo genérico (m, n) en la membrana circular tendremos m líneas nodales radiales y n círculos nodales, 
incluyendo el borde de la membrana. La vibración real será una superposición de los distintos modos. 
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2.4. Ondas estacionarias de presión en una caja rectangular 


Analicemos ahora las ondas estacionarias que se forman en un fluido contenido en una caja de paredes rectangulares de 
dimensiones a Xx b x c. Para ello obtendremos en primer lugar soluciones de la ecuación (1-73) mediante un método de 
separación de variables similar al utilizado en la sección 1.5.2, con la solución escrita como 


plz, y, 2,t) = X(DY(yYZ(IT() (2-15) 
que da una dependencia temporal similar a la ya obtenida para la ecuación (1-26), 
T(E) = AO 00D (1-27) 


y soluciones para la parte espacial similares a las de la ecuación (1-28) pero con una constante de separación más (k,), 
de forma que la solución es 


Ae itertky+k: 2+wt—a) > pe Ertut-a) 


plz, y, 2,t) = (2-16) 
con k = (ka, Ky, kz) el vector de onda y F'= (x, y, 2), siendo la relación entre las distintas constantes 
2 2 2 2 w? 
ka + hy + ko == (2-17) 


Si recordamos que el desplazamiento Y está relacionado con la presión mediante la ecuación (1-72) una onda como la 
anterior, con dependencia temporal en exp(¿wt) da, después de una doble integración en £, 


1 
pow? 


ú= Ve (2-18) 


ya que no aparecen términos convectivos (y así dd /dt = 04/0t). Como el desplazamiento debe anularse en las paredes 
de la caja, una solución de la forma! 


plz, y, 2,t) = pocos kx cos ky y cos TA A (2-19) 
nos daría unas componentes del desplazamiento 
Ug = -— cebo sen kx COS ky y COS ke 04-90) (2-20a) 
Uy = — e COS kh 1 sen ky y COS PO A (2-20b) 
Us == E cos kx cos ky y sen O (2-20c) 


que se hacen O para z = 0, y = 0 y z =0 y que también se anularán para x = a, y =b y z =Csi ky, k, y k, son de la 
forma 


kytai= TL 10,1, Ds: (-2la) 
kyb=1m m=0,1,2,... / (Q-21b) 
k0=xnn n=0,1, BD. (Q21c) 


que nos da unos valores para las frecuencias angulares de cada uno de los modos 


Ml e my? ny? 
w(l,m,n) = rc (2) + (5) +(2) ; (2-22) 


Obviamente al menos una de las l, m, n debe ser distinta de O. En este caso tendremos ! — 1 planos nodales perpendiculares 
al eje x en los que la componente x del desplazamiento es 0, m— 1 planos nodales perpendiculares al eje y y n— 1 planos 
nodales perpendiculares al eje z, que dividen el recinto para cada modo de vibración en l x m x n elementos de volumen 
prismáticos. La vibración más general será superposición lineal de los distintos modos, en el rango de validez del principio 
de superposición. 


TAI desarrollar la ecuación (2-15) aparecen otros sumandos que no satisfarán las condiciones de contorno, pero que también son solución de la 
ecuación de ondas sin condiciones en los límites. 
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